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Summary. Nuclear reactions transform atomic nuclei inside stars. This is the pro- 
cess of stellar nucleosynthesis. The basic concepts of determining nuclear reaction 
rates inside stars are reviewed. How stars manage to burn their fuel so slowly most 
of the time are also considered. Stellar thermonuclear reactions involving protons 
in hydrostatic burning are discussed first. Then I discuss triple alpha reactions in 
the helium burning stage. Carbon and oxygen survive in red giant stars because 
of the nuclear structure of oxygen and neon. Further nuclear burning of carbon, 
neon, oxygen and silicon in quiescent conditions are discussed next. In the subse- 
quent core-collapse phase, neutronization due to electron capture from the top of 
the Fermi sea in a degenerate core takes place. The expected signal of neutrinos 
from a nearby supernova is calculated. The supernova often explodes inside a dense 
circumstellar medium, which is established due to the progenitor star losing its out- 
ermost envelope in a stellar wind or mass transfer in a binary system. The nature 
of the circumstellar medium and the ejecta of the supernova and their dynamics are 
revealed by observations in the optical, IR, radio, and X-ray bands, and I discuss 
some of these observations and their interpretations. 

Keywords: nuclear reactions, nucleosynthesis, abundances; stars: interiors; 
supernovae: general; neutrinos; circumstellar matter; X-rays: stars 

1 Introduction 

The sun is not commonly considered a star and few would think of stars as 
nuclear reactors. Yet, that is the way it is, and even our own world is made 
out of the "fall-out" from stars that blew up and spewed radioactive debris 
into the nascent solar systerrn. 



Nuclear Astrophysics is the field concerning "the synthesis and Evolution 
of atomic nuclei, by thermonuclear reactions, from the Big Bang to the present. 

1 Lecture Notes on Kodai School on Synthesis of Elements in Stars; eds. Aruna 
Goswami & Eswar Reddy, Springer Verlag, 2009 
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What is the origin of the matter of which we are made?" pQ . Our high entropy 
universe resulting from the Big Bang, contains many more photons per particle 
of matter with mass, e.g. protons and neutrons. Because of the high entropy 
as the universe expanded, there was time to manufacture elements only upto 
helium and the major products of cosmic nucleosynthesis remained the light 
elements hydrogen and heliuuQ. Stars formed from this primordial matter. 
They used these elements, hydrogen and helium as fuel to generate energy like 
a giant nuclear reactojfl In the process, the stars could shine and manufacture 
higher atomic number elements like carbon, oxygen, calcium and iron which 
we and our world are made of. The heavy elements are either dredged up from 
the core of the star to the surface of the star from which they are dispersed by 
stellar wind or directly ejected into the interstellar medium when a (massive) 
star explodes. The Stardust is the source of heavy elements for new generation 
of stars and sun-like systems. 

Our sun is not a massive star. It burns hydrogen in a set of nuclear re- 
actions called the pp-chain, whereas the more massive stars presently burn 
their hydrogen by the so-called CNO-cycle0. Nevertheless, to put nuclear re- 
actions in stars in perspective, we shall start with a discussion of how these 
reactions proceed within the sun. There is a correspondence between the evo- 
lutionary state of a star, its external appearanc^f] and internal core conditions 
and the nuclear fuel it burns, - a sort of a mapping between the astronomers 
Hcrtzsprung-Russel diagram and the nuclear physicist's chart of the nuclide 
[4], until nuclear burning takes place on too rapid a time scale. 

The problem of evolution of stars to their explosion and subsequent in- 
teraction with the circumstellar medium has many time scales (ranging from 
tens of milliseconds to tens of thousands of years) and macroscopic length 
scales (from dimensions effectively that of a white dwarf to that of a super- 
nova remnant, i.e from few thousand kilometers to many tens of light years). 
The physics of supernova explosions is complex and diverse and in many parts, 
the explosion mechanism is still an unsolved problem. Even the constraining 
parameters and ingredients which makes the SN explode are still controver- 



2 Note however suggestions [5] [6] that early generation of stars called Pop III objects 
can also contribute to the abundance of 4 He seen in the universe today and the 
entire helium may not be a product of big bang nucleosynthesis alone. 

3 Our sun is slowly burning hydrogen into helium and is presently not exactly 
the same when it just started burning hydrogen. It will appear different once it 
exhausts all hydrogen it can burn in its core. In other words, nuclear reactions in 
stellar interiors determine the life-cycle of stars, apart from providing them with 
internal power for heat and light and manufacturing all the heavier elements that 
the early universe could not. 

4 Note however that CN cycle may have driven an early stage convection in the 
young Sun. 

Astronomers classify stars according to their colors and (absorption) line spectra 
and luminosities. Meghnad Saha showed [2j[3] the link between the classification 
scheme and temperature (and thermal ionization) of stellar atmosphere. 
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sial (see e.g. the discussion in [96 ). It is possible that the identification of 
the key aspects in the explosion may require seminal observations about the 
conditions in the supernova core other than the indirect evidence such as ex- 
plosion asymmetries, pulsar kicks or nucleosynthetic yields. Such observations 
may involve the future detection of strong neutrino signals and gravitational 
waves from a galactic supernova in future. Detectable neutrino signals from 
a supernova was seen in the case of SN 1987A, but since that target was in 
a neighboring satellite galaxy (the Large Magellanic Cloud), the number of 
neutrinos were too small to distinguish the fine points of theoretical issues. 

Since nuclear astrophysics is not usually taught at the master's level nu- 
clear physics specialization in our universities, these lecture notes are meant 
to be an introduction to the subject and a pointer to literature and Inter- 
net resource^. The emphasis in the first part of these lecture notes is on the 
nuclear reactions in the stars and how these are calculated, rather than how 
stars evolve. The latter usually forms a core area of stellar astrophysics. 

This article is organized essentially in the same sequence that a massive 
star burns successively higher atomic number elements in its core, until it col- 
lapses and explodes in a supernova. The introductory part discusses how the 
rates of thermonuclear reactions in (massive) stars are calculated, what the 
different classes of reactions are and how the stars (usually) manage to burn 
their fuels so slowl}0. The middle part describes the nuclear physics during 
the collapse phase of the massive star. The last part describes a few typical 
examples of what can be learned by optical, IR and X-ray studies about nucle- 
osynthesis and dynamics of explosion in supernovae and supernova remnants 
such as Cassiopeia A, SN 1987A etc. Only core-collapse supernovae are dis- 
cussed in these lectures, those that arise from massive stars (e.g. stars more 
massive than 8M with typical solar metallicity at the time they start burn- 

6 See for example, [7] for a course of nuclear astrophysics, and the International 
Conference Proceedings under the title: "Nuclei in the Cosmos" for periodic re- 
search conferences in the field. Valuable nuclear astrophysics datasets in machine 
readable formats can be found at sites: [8], [9|. A new and updated version of the 
nuclear reactions rate library REACLIB for astrophysics is now being maintained 
as a public, web-based version at [10] . A complementary effort to develop soft- 
ware tools to streamline the process of getting the latest and best information 
into this new library is available at [11] (see \1'2\). Much of the material dis- 
cussed in the first part of these notes can be found in textbooks in the subject, 
see e.g. [13] . [14] . [4], [15], [1] etc. There is also a recent book on the subject 
by Richard Boyd [TB] that among other topics describes terrestrial and space 
born instruments operating in service to nuclear astrophysics. A Workshop on 
Solar fusion cross sections for the pp chain and CNO cycle held in 2009 by 
the Institute of Nuclear Theory is expected to result in a Reviews of Modern 
Physics article. Supernovae of various types are the sites where nuclear reactions 
in stars or explosions are of prime importance. For Internet resources to two 



recent Schools on these topics, see 


http://icts.tifr.res.in/sites/Sgrb/Programme 


and http: / / www.tifr.res.in/~sn2004 





7 These issues were discussed in an earlier SERC School |17] . 
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ing hydrogen in their cores, i.e. at the "Zero Age Main Sequence" , before any 
mass was lost from their surface). We shall not discuss the type la SN^| in 
these lectures. Abundance of elements in our galaxy Milky Way give impor- 
tant information about how stars affect the element and isotopic evolution in 
various parts of the galaxy. For a study of the evolution of elements from C 
to Zn in the galactic halo and disk, with yields of massive stars and type la 
SNe, see [20] and elsewhere in these Proceedings. 

2 Stars and their thermonuclear reactions 

While referring to Sir Ernest Rutherford "breaking down the atoms of oxy- 
gen and nitrogen" , Arthur Eddington remarked: "what is possible in the 
Cavendish Laboratory may not be too difficult in the sun" [21]. Indeed this 
is the case, but for the fact that a star does this by fusion reactions, rather 
than by transfer reactions, - in the process giving out heat and light and 
manufacturing fresh elements. Of all the light elements, hydrogen is the most 
important one in this regard, because: a) it has a large universal abundance, 
b) considerable energy evolution is possible with it because of the large bind- 
ing energies of nuclei that can be generated from its burning and c) its small 
charge and mass allows it to penetrate easily through the potential barriers 
of other nuclei. A star burns its fuel in thermonuclear reactions in the core 
where the confinement of the fuel is achieved in the star's own gravitational 
field. These reactions remain "controlled" , or self-regulatecd, as long as the 
stellar material remains non-degenerate. 

8 The supernovae are classified by astronomers on the basis of their optical spectra 
at the time of their maximum light output. Those that do not show the presence 
of hydrogen in their spectra are classified as type I SNe. A subclass of them, type 
la's are believed to arise from thermonuclear explosions in the electron degenerate 
cores of stars less massive than 8Mq and are very useful to map the geometry of 
our universe, because they serve as calibratable "standard" candles. These "ther- 
monuclear supernovae" are usually more luminous in the optical bands than the 
core-collapse varieties, but while the former class put out several MeVs of energy 
per nucleon, the core-collapse SNe or ccSNe, emit several hundreds of MeVs per 
nucleon from their central engines (mostly in down-scattered neutrinos). Apart 
from the missing hydrogen lines, the type la SNe show an absorption "trough" 
in their spectra around maximum light at 6150 A, due to blue shifted Si II lines 
[18] . Of the other type I SNe which show no hydrogen and no Si II trough, some 
show helium in their early spectra (type lb) and others have little or no helium 
and have a composition enhanced in oxygen (type Ic) |19| . These, (lb and Ic) 
together with the type lis constitute the core collapse SNe. 

9 There are however examples to the contrary when thermonuclear reactions take 
place in an explosive manner, e.g. when a whole white dwarf (resulting from an 
evolved intermediate mass star) undergoes merger with another and explodes, as 
nuclear fuel (carbon) is ignited under degenerate conditions, such as in a type la 
supernova; explosive thermonuclear reactions also take place in post-bounce core 
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The recognition of the quantum mechanical tunneling effect prompted 
Atkinson and Houtermans [55] to work out the qualitative treatment of en- 
ergy production in stars. They suggested that the nucleus serves as both a 
cooking pot and a trap. Binding energy difference of four protons and two 
electrons (the nuclear fuel) and their ash, the helium nucleus, some 26.7 MeV 
is converted to heat and light that we receive from the surF°l. The photons in 
the interior are scattered many a times, for tens of millions of years, before 
they reach the surface of the sun. Some of the reactions produce neutrinos, 
which because of their small cross-section for interaction, are not stopped or 
scattered by overlying matter, - but stream out straight from the core. Neu- 
trinos are thus the best probes of the stellar core [23] 1106] , while the photons 
bear information from their surface of last scattering - the photosphere. 

2.1 Why do the stars burn slowly: a look at Gamow peaks 

The sun has been burning for at least 4.6 billion year^|. How does it manage 
to burn so slowl£3? Under the ambient conditions in the core, the relevant 



collapse supernovae, when the hydrodynamic shock ploughs through unburnt Si- 
or O-layers in the mantle. 

10 Lord Kelvin surmised in the nineteenth century that the solar luminosity is 
supplied by the gravitational contraction of the sun. Given the solar luminosity, 
this immediately defined a solar lifetime (the so-called Kelvin-Helmholtz time): 
tkh = GMq/ RqLq ~ few x 10 7 yr. This turned out to be much shorter than 
the estimated age of the earth at that time known from fossil records, and led to 
a famous debate between Lord Kelvin and the fossil geologist Cuvier. In fact, as 
noted above modern estimates of earth's age are much longer and therefore the 
need to maintain sunshine for such a long time requires that the amount of energy 
radiated by the sun during its lifetime is much larger than its gravitational energy 
or its internal (thermal) energy: Lq x tu[ c S> GMq/ Rq. This puzzle was resolved 
only later with the realization that the star can tap its much larger nuclear energy 
reservoir in its core through thermonuclear reactions. The luminosity of the sun 
however is determined by an interplay of atomic and gravitational physics that 
controls the opacity, chemical composition, the balance of pressure forces against 
gravity, etc. Nuclear physics determines how fast nuclear reactions go under under 
feedback control determined by the ambient conditions. 

11 Lord Rutherford [24] determined the age of a sample of pitchblende, to be 700 mil- 
lion years, by measuring the amount of uranium and radium and helium retained 
in the rock and by calculating the annual output of alpha particles. The oldest 
rock found is from Southwest Greenland: « 3.8 Gyr old [15| . Radioactive dating 
of meteorites point to their formation and the solidification of the earth about 
4.55 ±0.07 years ago [25| . Since the sun and the solar system formed only slightly 
before, their age at isolation and condensation from the interstellar medium is 
taken to be 4.6 Gyr [25] . 

12 The Nobel prize citation of Hans Bethe (1967) who solved this problem, noted 
that this " concerns an old riddle. How has it been possible for the sun to emit 
light and heat without exhausting its source not only during the thousands of 
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thermonuclear reaction cross sections are very small^j. For reactions involving 
charged particles, nuclear physicists often encounter cross-sections near the 
Coulomb barrier of the order of millibarns. One can obtain a characteristic 
luminosity Lc based on this cross section and the nuclear energy released per 
reaction [15] : 

L c ~ eNAE/rc (1) 

where e w 10 -2 is the fraction of total number of solar nuclei N ~ 10 57 
that take part in nuclear fusion reactions generating typically AE ~ 25 MeV 
in hydrogen to helium conversion. Here, the tq is the characteristic time scale 
for reactions, which becomes minuscule for the cross-sections at the Coulomb 
barrier and the ambient density and relative speed of the reactants etc: 

1 lCT 8 s 

" nav [n/(10 26 cm~ 3 )][er/l mbarn][v/10 9 cm s _1 ] 

This would imply a characteristic luminosity of L c 10 20 Lq, even for a 
small fraction of the solar material taking part in the reactions (i.e. e ~ 10~ 2 ). 
If this was really the appropriate cross-section for the reaction, the sun would 
have been gone very quickly indeed. Instead the cross-sections are much less 
than that at the Coulomb barrier penetration energy (say at proton energies 
of 1 MeV), to allow for a long lifetime of the sun (in addition, weak-interaction 
process gives a smaller cross-section for some reactions than electromagnetic 
process, - see Section 3.1). 

Stellar nuclear reactions can be either: a) charged particle reactions (both 
target and projectile are nuclei) or b) neutral particle (neutron) induced re- 
actions. Both sets of reactions can go through either a resonant state of an 
intermediate nucleus or can be a non-resonant reaction. In the former reac- 
tion, the intermediate state could be a narrow unstable state, which decays 
into other particles or nuclei. In general, a given reaction can involve both 



centuries the human race has existed but also during the enormously long time 
when living beings needing the sun for their nourishment have developed and 
flourished on our earth thanks to this source? The solution of this problem seemed 
even more hopeless when better knowledge of the age of the earth was gained. 
None of the energy sources known of old could come under consideration. A very 
important part of his work resulted in eliminating a great number of thinkable 
nuclear processes under the conditions at the center of the sun, after which only 

two possible processes remained (Bethe) attempted a thorough analysis of 

these and other thinkable processes necessary to make it reasonably certain that 
these processes, and only these, are responsible for the energy generation in the 
sun and similar stars. " 

This makes the experimental verification of the reaction cross-sections a very 
challenging task, requiring in some cases, extremely high purity targets and pro- 
jectiles so that relevant small event rates are not swamped by other reaction 
channels and products (see Rolfs and Rodney, Chapter 5 13 ). 



Thermonuclear fusion and collapse and explosion in massive stars 7 

types of reaction channels. In charged particle induced reactions, the cross- 
section for both reaction mechanisms drops rapidly with decreasing energy, 
due to the effect of the Coulomb barrier (and thus it becomes more difficult 
to measure stellar reaction cross-sections accurately). In contrast, the neu- 
tron induced reaction cross-section is very large and increases with decreasing 
energy (here, resonances may be superposed on a smooth non-resonant yield 
which follows the l/v ~ \j\f~E dependence). These reaction rates and cross- 
sections can be then directly measured at stellar energies that are relevant (if 
such nuclei are long lived or can be generated). The schematic dependence 
of the cross-sections are shown in Fig. [TJ We shall not discuss the neutron 
capture elements in these notes. 




INTERACTION ENERGY E 



Fig. 1. Dependence of total cross-sections on the interaction energy for neutrons 
(top panel) and charged particles (bottom panel). Note the presence of resonances 
(narrow or broad) superimposed on a slowly varying non-resonant cross-section (af- 
ter [13]). 



2.2 Gamow peak and the astrophysical S-factor 

The sun and other "main-sequence" stars (burning hydrogen in their core qui- 
escently) evolve very slowly by adjusting their central temperature such that 
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the average thermal energy of a nucleus is small compared to the Coulomb 
repulsion an ion-ion pair encounters. This is how stars can live long for astro- 
nomically long times. A central temperature T > 10 7 K (or T 7 > 1, hereafter 
a subscript x to a quantity, indicates that quantity in units of 1CP) is required 
for sufficient kinetic energy of the reactants to overcome the Coulomb barrier 
and for thermonuclear reactions involving hydrogen to proceed at an effective 
rate, even though fusion reactions have positive Q values i.e. net energy is lib- 
erated out of the reactions. The classical turning point radius for a projectile 
of charge Z 2 and kinetic energy E p (in a Coulomb potential Vc — ZiZ 2 e 2 /r, 
and effective height of the Coulomb barrier Ec = Z\Z 2 e? j R n = 550 keV for a 
p + p reaction), is: r c i — ZiZ 2 e 2 /E p . Thus, classically a p -I- p reaction would 
proceed only when the kinetic energy exceeds 550 keV. Since the number of 
particles traveling at a given speed is given by the Maxwell Boltzmann (MB) 
distribution (f>(E), only the tail of the MB distribution above 550 keV is effec- 
tive when the typical thermal energy is 0.86 keV ( Tg = 0.01). The ratio of the 
tails of the MB distributions: 0(550 keV)/0(O.86 keV) is quite minuscule, and 
thus classically at typical stellar temperatures this reaction will be virtually 
absent. 

Although classically a particle with projectile energy E p cannot penetrate 
beyond the classical turning point, quantum mechanically, one has a finite 
value of the squared wave function at the nuclear radius R n : |^>(i2„)| 2 . The 
probability that the incoming particle penetrates the barrier is: 

WR C )\ 2 (3) 

where ip( r ) are the wave- functions at corresponding points. Bethe |27j 
solved the Schroedinger equation for the Coulomb potential and obtained the 
transmission probability:- 

P ( OT^R r ^HRc/Rn-l) 172 R n] \ 

with K = [2fi/h 2 (E c — E)] 1 / 2 . This probability reduces to a much simpler 
relation at the low energy limit: E ^ E c , which is equivalent to the clas- 
sical turning point R c being much larger than the nuclear radius R n . The 
probability is: 



P = exp(-27rr/) = exp[-27rZiZ 2 e 2 /(hv)] = exp[-31.3Z 1 Z 2 (^) 1 / 2 ] (5) 

E 

where in the second equality, \i is the reduced mass in Atomic Mass Units 
and E is the center of mass energy in keV. The exponential quantity involving 
the square brackets in the second expression is called the "Gamow factor". 
The reaction cross-section between particles of charge Z\ and Z 2 has this 
exponential dependence due to the Gamow factor. In addition, because the 
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cross-sections are essentially "areas": proportional to 7r(A/27r?i) 2 oc l/E, it 
is customary to write the cross-section, with these two energy dependences 
filtered out: 

«{E) = C -^^S(E) (6) 

where the factor S(E) is called the astrophysical S-factor. The S-factor 
may contain degeneracy factors due to spin, e.g. [(2 J+ l)/(2Ji + 1)(2 J 2 + 1)] 
as reaction cross-sections are summed over final states and averaged over ini- 
tial states. Because the rapidly varying parts of the cross-section (with energy) 
are thus filtered out, the S-factor is a slowly varying function of center of mass 
energy, at least for the non-resonant reactions. It is thus much safer to ex- 
trapolate S(E) to the energies relevant for astrophysical environments from 
the laboratory data, which is usually generated at higher energies (due to 
difficulties of measuring small cross-sections) , than directly extrapolating the 
cr(E), which contains the Gamow transmission factor (see Fig. [5]). Addition- 
ally, in order to relate a(E) and 3(E), quantities measured in the laboratory 
to these relevant quantities in the solar interior, a correction factor /o due to 
the effects of electron screening needs to be taken into account [28] . 

In the stellar core with a temperature T, reacting particles have many 
different velocities (energies) according to a Maxwell Boltzmann distribution:- 

cx E 1/2 exp[-E/kT] (7) 

Nuclear cross-section or the reaction rates which also depend upon the 
relative velocity (or equivalently the center of mass energy) therefore need to 
be averaged over the thermal velocity (energy) distribution. Therefore, the 
thermally averaged reaction rate per particle pair is: 

f°° 8 1 /*°° 

< av >= / ct>(v)a(v)vdv = (— ) l/ —— / a(E)E exp(-E/kT)dE 
Jo ) ' Jo 

(8) 

The thermally averaged reaction rate per pair is, utilizing the astrophysical 
S-factor and the energy dependence of the Gamow-factor: 

< " >= o ,,2 ww r s<E,exp[ - s - < 9) 

with b 2 = Eq = 2^(vre 2 Z^/fi) 2 = 0.978/iZfZ 2 MeV, E G being called 
the Gamow energy. Note that in the expression for the reaction rate above, 
at low energies, the exponential term exp(— b/\/E) = exp(— y^Ec/E)) be- 
comes very small whereas at high energies the Maxwell Boltzmann factor 
E 1 / 2 exp(— E/kT) vanishes. Hence there would be a peak (at energy, say, Eq) 
of the integrand for the thermally averaged reaction rate per pair (see Fig. [3j. 
The exponential part of the energy integrand can be approximated as: 



(v) = iwv" 



2nkT , 







J exp 


~ 2kT 
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Energy E — — 



Fig. 2. Cross-section and astrophysical S-factor for charged particle reactions as a 
function of beam energy. The effective range of energy in stellar interiors is usually 
far less than the Coulomb barrier energy Ec or the lower limit El of laboratory 
measurements. The y-scale is logarithmic for cross-section but linear for S-factor; 
thus the cross section drops sharply in regions of astrophysical interest, while the 
change is much less severe for the S-factor. The extrapolation of laboratory data to 
lower energies relevant for astrophysical situations is more reliable for S-factor. 



cxp[- j^-bE- 1 * 72 ] -Cexp 



,E-Eov 2 ' 
' A/2 ' 



(10) 



where 



C = cxp(-E /kT - bE 1/2 ) = cxp(-3E /kT) = exp(-r) 
E = (bkT/2)i = 1.22keV(ZfZ^T|)3 



A = 4(E kT/3)i = 0.75keV(Z 2 Z 2 ATg)5 



Since most stellar reactions happen in a fairly narrow band, of energies, 
S(E) will have a nearly constant value over this band averaging to So- With 
this, the reaction rate per pair of particles, turns out to be: 
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E 

Energy E — ^- 

Fig. 3. The Gamow peak is a convolution of the energy distribution of the Maxwell 
Boltzmann probability and the quantum mechanical Coulomb barrier transmission 
probability. The peak in the shaded region near energy Eq is the Gamow peak that 
gives the highest probability for charged particle reactions to take place. Usually the 
Gamow peak is at a much higher energy than kT, and in the figure the ordinate scale 
(for the Gamow peak) is magnified with respect to those of the M-B and barrier 
penetration factors. See also Table [T] 



<<7V>= [All^ [°° e -r- ^(^ a ) 2 ^ = 4 .51Q 14 - r |Vr 2 e- T cm 3 s- 1 

L fcT J TT/J, J AZiZ 2 

(11) 

Here, 

T = 3E /kT = 42.5(Z%Z%ft/T 6 )* (12) 
The maximum value of the integrand in the above equation is: 

Imax = exp(-r) 

The values of Eo,I max , A, etc., apart from the Coulomb barrier for several 
reactions are tabulated in Table Q] for T 6 = 15. 

As the nuclear charge increases, the Coulomb barrier increases, and the 
Gamow peak Eq also shifts towards higher energies. Note how rapidly the 
maximum of the integrand I max decreases with the nuclear charge and the 
Coulomb barriers. The effective width A is a geometric mean of Eo and kT, 
and A/2 is much less rapidly varying between reactions (for kT <C Eo). The 
rapid variation of I max indicates that of several nuclei present in the stel- 
lar core, those nuclear pairs will have the largest reaction rates, which have 
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the smallest Coulomb barrier. The relevant nuclei will be consumed most 
rapidly at that stage. (Note however that for the p+p reaction, apart from 
the Coulomb barrier, the strength of the weak force, which transforms a pro- 
ton to a neutron also comes into play) . 

When nuclei of the smallest Coulomb barrier are consumed, there is a 
temporary dip in the nuclear generation rate, and the star contracts gravi- 
tationally until the temperature rises to a point where nuclei with the next 
lowest Coulomb barrier will start burning. At that stage, further contraction is 
halted. The star goes through well defined stages of different nuclear burning 
phases in its core dictated by the height of the Coulomb barriers of the fuels. 
Note also from the Table [TJ how far Eq, the effective mean energy of reaction 
is below the Coulomb barrier at the relevant temperature. The stellar burning 
is so slow because the reactions are taking place at such a far sub-Coulomb 
region, and this is why the stars can last so long. 

Table 1. Parameters of the thermally averaged reaction rates at Te = 15. 



Reaction 


Coulomb Gamow 


I max 


A 






Barrier 


Peak (Eo) (e- SE ^ kT ) 








(MeV) 


(keV) 




(keV) 




P + P 


0.55 


5.9 


1.1 x 10~ 6 


6.4 


7 x 10" 6 


P + N 


2.27 


26.5 


1.8 x 10~ 27 


13.6 


2.5 x 10 -26 


a + C 12 


3.43 


56 


3 x 10~ 57 


19.4 


5.9 x 10~ 56 


Ql6 + Ql6 


14.07 


237 


6.2 x 10" 239 


40.4 


2.5 x 10" 237 



The above discussion assumes that a bare nuclear Coulomb potential is 
seen by the charged projectile. For nuclear reactions measured in the labora- 
tory, the target nuclei are in the form of atoms with electron cloud surrounding 
the nucleus and giving rise to a screened potential - the total potential then 
goes to zero outside the atomic radius. The effect of the screening is to reduce 
the effective height of the Coulomb barrier. Atoms in the stellar interiors are 
in most cases in highly stripped state, and nuclei are immersed in a sea of 
free electrons which tend to cluster near the nucleus. When the stellar den- 
sity increases, the so called Debye-Huckel radius Rd = (kT/Ane 2 pA^f) 1 / 2 , 
(here: £ = J2i(^f + Zi)Xi/Ai) which is a measure of this cluster "radius", 
decreases, and the effect of shielding upon the reaction cross-section becomes 
more important. This shielding effect enhances thermonuclear reactions in- 
side the star. The enhancement factor fo = exp(0.188ZiZ2£p 1 / 2 T 6 3 ^ 2 , varies 
between 1 and 2 for typical densities and compositions [5H] but can be large 
at high densities. 
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3 Hydrogen burning: the pp chain 

The quantitative aspects of the problem of solar energy production with de- 
tails of known nuclear physics of converting hydrogen into helium was first 
worked out by von Weizsacker (1937-38) [32], [33] and Bethe & Critchfield 
(1938-1939) [34], which made it clear that two different sets of reactions : the 
p-p chains and the CN cycle can do this conversion. This happens in the core 
of the star initially (at the "main sequence" stage), and then later in the life 
of a star in a shell of burning hydrogen around an inert core of He. 

In the first generation of stars in the galaxy only the p-p cycle may have 
operated. In second generation, heavier elements like C, N from the ashes 
of burning in previous stars are available and they too can act as catalysts 
to have thermonuclear fusion of hydrogen to helium. Since in the very first 
generation, the heavier elements like C, N were practically absent, all stars 
including the massive stars, burnt hydrogen through the p-p cycle. [A recent 
discovery ( 35J) of a low-mass star with an iron abundance as low as 1/200,000 
of the solar value (compare the previous record of lowest iron abundance less 
than 1/10,000 that of the sun), suggests that such first generation stars are 
still around]. 

The sun with a central temperature of 15.7 million degrees, (Tg = 15.7) 
burns by p-p chains. Slightly more massive star (with central temperature 
2e > 20) burns H by the CNO cycle also. Davis et al.s' solar neutrino exper- 
iment [23], which in 1968 had only an upper limit of the neutrino flux, itself 
put a limit of less than 9% of the sun's energy is produced by the carbon- 
nitrogen cycle (the more recent upper limit [36j is 7.3%, from an analysis of 
several solar neutrino experiments, including the Kamland measurements). 
Note however that for the standard solar model, the actual contribution of 
CNO cycle to solar luminosity is ~ 1.5% [IS]). In CNO cycle, nuclei such as 
C, N, O serve as "catalysts" do in a chemical reaction. The pp-chain and the 
CNO cycle reaction sequences are illustrated in Figs. 141 and [Till 

The pp-chain begins with the reaction p + p^d + c + + ^ c . Bethe and 
Critchfield [34] showed that weak nuclear reaction is capable of converting 
a proton into a neutron during the brief encounter of a scattering event. 
(This reaction overcomes the impasse posed by the instability of 2 He in the 
p + p — > 2 He reaction via the strong or electromagnetic interactions, and as 
the next nucleus 3 Li reachable via these interactions is also unstable as a 
final product). Since a hydrogen atom is less massive than a neutron, such a 
conversion would be endothermic (requiring energy), except for the fact that 
a neutron in a deuterium nucleus 2 D can form a bound state with the proton 
with a binding energy of 2.224 MeV - thus making the reaction exothermic 
with an available kinetic energy of 0.42 MeV. The released positron soon pair 
annihilates into photons making the total energy released to be 1.442 MeV. 

Because of the low Coulomb barrier, in the p+p reaction (E c — 0.55 MeV), 
a star like the sun would have consumed all its hydrogen quickly (note the 
relatively large value of (A)I max in TableQ}, were it not slowed down by the 
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THE REACTIONS Of THE P-P CHAIN,4p— 4 He + 2e + + 2*,+!}^ (BP2000 model) 

( l y ■ 0.26 M«V 
p(p. «M d j Fi/ « exIO^cm"^- 1 
<J (P.7) 3 H« 
85% I 15% 

I I 

3 H»( 3 He,2p) 4 He 3 He(n.r) 7 Be 

1 5% | 0.0 2% 

I I 

7 Be(«-,»0 7 U 7 fe(p.7) 8 B , t _ ?3 ^ 
7 Li (p. «) 4 He 8b (#,*) 8 Be* j ^. 5.8x1 0+ 8 cm" 2 s" 1 
V («) *He 

CHAIN I CHAIN II CHAIN 1 

Qgff « 26.20 MaV Q eff - 25.66 M«V Q eff - 19.17 MeV 
( 2.0% loss) ( 4.0% loss) ( 28.3% loss) 

Fig. 4. The p-p chain starts with the formation of deuterium and 3 He. Thereafter, 
3 He is consumed in the sun 85% of the time through ppl chain, whereas ppll and 
ppIII chains together account for 15% of the time in the Bahcall Pinsonneault 2000 
solar model. The ppIII chain occurs only 0.02% of the time, but the 8 B /3 + -decay 
provides the higher energy neutrinos (average E v — 7.3 MeV). The net result of 
the chains is the conversion of four protons to a helium, with the effective Q- values 
(reduced from 26.73 MeV) as shown, due to loss of energy in escaping neutrinos. See 
|38II37] for updated branching ratios and neutrino fluxes for BPS2008(AGS) model. 

weakness of the weak interactions. The calculation of probability of deuteron 
formation consists of two separate considerations: 1) the penetration of a 
mutual potential barrier in a collision of two protons in a thermal bath and 2) 
the probability of the /3-decay and positron and neutrino emission. Bcthe and 
Critchfrcld used the original Fermi theory (point interaction) for the second 
part, which is adequate for the low energy process. 

3.1 Cross-section for deuterium formation 

The total Hamiltonian H for the p-p interaction can be written as a sum of 
nuclear term H n and a weak- interaction term H w . As the weak interaction 
term is small compared to the nuclear term, first order perturbation theory 
can be applied and Fermi's "Golden rule" , gives the differential cross-section 
as: 



E„- 0.81 MeV 

F„« 0.47x1 10 cm~V 



Thermonuclear fusion and collapse and explosion in massive stars 



15 
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here p(E) — dN/dE, is the density of final states in the interval dE and 
Vi is the relative velocity of the incoming particles. For a given volume V, the 
number of states dn between p and p+dp is:- 

dN = dn e dn v = (V^^)(V^pL) (14) 

By neglecting the recoil energy of deuterium (since this is much heavier 
than the outgoing positron in the final state) and neglecting the mass of the 
electron neutrino, we have: 

E = E e + E v = E e + cp v (15) 
and dE = dE v — cp v , for a given E e and, 



16n 2 V 2 

P (E) = dN(E)/dE = dn e {dn v /dE) = p 2 e (E - E e f dp e = p(E e )dp e 

(16) 

The matrix element that appears in the differential cross section, may be 
written in terms of the initial state wave function W% of the two protons in the 
entrance channel and the final state wave function <?y as: 

H if = J [V&VrfHpWidT (17) 

If the energy of the electron is large compared to Z x Rydberg (Rydberg 
Roc = 2ir 2 me /c/i 3 ), then a plane wave approximation is a good one: = 
1/ (v / V")cxp(ifc e .r) where the wave-function is normalized over volume V. (For 
lower energies, typically 200 keV or less, the electron wave-function could be 
strongly affected by nuclear charge (see [39])). Apart from this, the final state 
wave function: [tf^^e^i/] has a deuteron part ^ whose radial part rapidly 
vanishes outside the nuclear domain (Rq), so that the integration need not 
extend much beyond r ~ Rq (for example, the deuteron radius R4 = 1.7 
fm). Note that because of the Q- value of 0.42 MeV for the reaction, the 
kinetic energy of the electron (K e < 0.42 MeV) and the average energy of 
the neutrinos [E v = 0.26 MeV) are low enough so that for both electrons and 
neutrino wave-functions, the product kRo < 2.2 x 10~ 3 and the exponential 
can be approximated by the first term of the Taylor expansion: 

# e = l/(W)[l + i(fee.r)] ~ 1/(VV) (18) 

and 

Then the expectation value of the Hamiltonian, for a coupling constant g is: 
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Hif = J [*d*c*,]*H^*idr = 1 J [* d ]**idr (19) 
The dr integration can be broken into space and spin parts M space and 

2n 167r 2 

da = —-^g 2 M 2 spin M% ac j>l{E - E e fd Pe (20) 
The total cross-section upto an electron energy of E is proportional to: 



E (™ ^5 r W 

p 2 e {E-E e) ,, r , - 

it c J i 

J- 



2^ ^f dpe= ^L {W ? - if'^W - W e ) 2 W e dW e (21) 



where W = (E + m e c ) / m e c ■ The integral over W can be shown as: 



f(W) = (W 2 - Ifl^ - ^ - 1] + ^ ln[W + (W 2 - I) 1 / 2 ] (22) 



so that: 

5 4 

mzc 



o Ml f{W)g'Mi pace Mi pin (23) 
At large energies, the factor f(W) behaves as: 

f(W) cx W 5 cx ^E 5 (24) 

The final state nucleus (deuterium in its ground state) in the reaction: 
p + p — > d + c + + ^ c , has = 1 + , with a predominant relative orbital angular 
momentum If — and Sf = 1 (triplet S-state). For a maximally probable 
super-allowed transition, there is no change in the orbital angular momentum 
between the initial and final states of the nuclei. Hence for such transitions, 
the initial state two protons in the p + p reaction must have h — 0. Since the 
protons are identical particles, Pauli principle requires Si = 0, so that the total 
wave-function will be antisymmetric in space and spin coordinates. Thus, we 
have a process: \Si = 0, k — >^ \Sf = 1, // = >. This is a pure Gamow- 
TelleiO transition with coupling constant g = Ca (the axial vector coupling 
component can be obtained, from the pure GT decay 6 He(0 + ) Li(l + )). 



14 In the beta-decay allowed approximation, we neglect the variation of the lepton 
wave-functions over the nuclear volume and the nuclear momentum (this is equiv- 
alent to neglecting all total lepton orbital angular momenta L > 0). The total an- 
gular momentum carried off by the leptons is their total spin: i.e. S = 1 or 0, since 
each lepton has When the lepton spins in the final state are anti-parallel, 
s E + s v = stot = the process is the Fermi transition with Vector coupling con- 
stant g = Cv (e.g. a pure Fermi decay: 14 0(Jf = 0+) -^ 14 N(JJ = 0+)). When 
the final state lepton spins are parallel, s? + s„ = s to t = 1, the process is Gamow- 
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The spin matrix element in the above expression for energy integrated 
cross-section a, is obtained from summing over the final states and averaging 
over the initial states and dividing by 2 to take into account that we have two 
identical particles in the initial state. Thus, 

1 777 5 r 4 M 2 M 2 ■ 

X = - = VJ^^ W ^ 2 ^ ( 25 ) 

t 2n 6 n Vi 2 

where, M 2 pin = @ j^i^h+i) = ^' ^ n< ^ s P acc matrix element is: 

/>OC 

M space = / X f{r) Xl {r)r 2 dr (26) 
Jo 

in units of cm 3 / 2 . The above integral contains the radial parts of the nu- 
clear wave- functions x(r), and involves Coulomb wave- functions for barrier 
penetration at (low) stellar energies. The integral has been evaluated by nu- 
merical methods ([42 ), and Fig. [5] shows schematically how the M space is 
evaluated for the overlap of the deuterium ground state wave-function with 
the initial pair of protons state. (See also [13], [S] for details of calculations of 
the overlap integral and writing the astrophysical S-factor in terms the beta 
decay rate of the neutron [43) which takes into account of radiative correc- 
tions to the axial- vector part of the neutron decay through an effective matrix 

Teller with g — Ca- For Fermi coupling, there is no change in the (total) angular 
momentum between the initial and final states of the nuclei {A J = 0). For the 
Gamow- Teller coupling, the selection rules are: A3 = or ±1 (but the possibility 
AJ = is excluded between two states of zero angular momentum). The size of 
the matrix element for a transition depends on the overlap of the wave-functions 
in the initial and final states. In the case of "mirror pair" of nuclei (the nucleus 
Az = (2Z + l)z is the mirror of the nucleus (2Z + l)z+i), the wave-functions 
are very much alike as shown through simple heuristic arguments [40] . For these 
nuclei, ft-vahies range from ~ 1000 — 5000 and are called super-allowed tran- 
sitions. For super-allowed transitions, which have maximum decay probabilities, 
there are no changes in the orbital angular momentum between the initial and 
final states of the nuclei. In the p + p^D+e + +f c reaction, the initial proton 
state is antisymmetric to an interchange of space and spin coordinates and the 
final deuteron is symmetric in this respect. (In fact when the two protons are 
in the S state (which is most favorable for their close approach), their spins will 
be anti-parallel (a singlet state) whereas the ground state of the deuteron is a 
triplet S state). If this were the complete description of the exchange symmetry 
properties of the Gamow- Teller transition (permitting a change of spin direction 
of the proton as it transforms to a neutron, changing the total spin by one unit) 
advocated here this would actually be forbidden. However in the use of configu- 
ration space in beta-decay process one must include isotopic spin as well. The S 
state of the two protons is symmetric to exchange of this coordinate, whereas the 
deuteron (consisting of both a proton and a neutron) function is antisymmetric 
in this coordinate. In the complete coordinate system the transition is from an 
initial antisymmetric state to another antisymmetric final state accompanied by 
a positron emission ([41] i. 
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element, the assumption being that these are the same as that for the proton 
beta decay in the pp reaction above). In the overlap integral one needs only 
the S-wave part for the wave- function of the deuteron ip^, as the D-wave part 
makes no contribution to the matrix element [42] , although its contribution 
to the normalization has to be accounted for. The wave-function of the initial 
two-proton system ip p is normalized to a plane wave of unit amplitude, and 
again only the S-wave part is needed. The asymptotic form of tp p (well outside 
the range of nuclear forces) is given in terms of regular and irregular Coulomb 
functions and has to be defined through quantities related to the S-wave phase 
shifts in p-p scattering data). The result is a minuscule total cross-section of 
a = 10 _47 cm 2 at a laboratory beam energy of E p — 1 MeV, which cannot be 
measured experimentally even with milliampere beam currents. 

The reaction p + p — > d + e + + v c is a non-resonant reaction and at all 
energies the rate varies smoothly with energy (and with stellar temperatures) , 
with 5(0) = 3.8 x l(r 22 keV barn and dS(Q)/dE = 4.2 x l(r 24 barn. At for 
example, the central temperature of the sun Tg = 15, this gives: < av > PP — 
1.2 x 10 -43 cm 3 s _1 . For density in the center of the sun p = 100 gm cm -3 and 
and equal mixture of hydrogen and helium (Xh = Xn e = 0.5), the mean life 
of a hydrogen nucleus against conversion to deuterium is th(H) = X/Nh < 
av > PP r ^ 10 10 yr. This is comparable to the age of the old stars. The reaction 
is so slow primarily because of weak interactions and to a lesser extent due 
to the smallness of the Coulomb barrier penetration factor (which contributes 
a factor ~ 10~ 2 in the rate), and is the primary reason why stars consume 
their nuclear fuel of hydrogen so slowly. For a calculation of the weak capture 
of protons on protons using calculated wavefunctions obtained from modern, 
realistic high precision interactions, see [29] 

3.2 Deuterium burning 

Once deuterium is produced in the weak interaction mediated p + p reaction, 
the main way this is burnt in the sun turns out to be: 

d + p^ 3 He +7 (27) 

This is a non-resonant direct capture reaction to the 3 He ground state with 
a Q-value of 5.497 MeV and S(Q) = 2.5 x 10~ 3 keV barn. The angle averaged 
cross-sections measured as a function of proton + deuterium center of mass 
energy, where the capture transitions were observed in gamma-ray detectors 
at several angles to the incident proton beam direction, are well explained by 
the direct capture model (see Fig. [6] after [13]). The LUNA collaboration [31] 
has measured the cross section down to a 2.5 keV cm. energy, well below the 
solar Gamow peak using a segmented Bismuth germanate (BGO) gamma-ray 
detector and found the S(E) factor to be in fair agreement with extrapolation 
of data at higher energies. 

The reactions comprising the rest of the (three) pp-chains start out with 
the predominant product of deuterium burning: 3 He (manufactured from d+p 
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Fig. 5. Schematic representation (after [13] ) of the numerical calculation of the 
spatial part of the matrix element M space in the p + p ^ d + e + + ^ c reaction. The 
top part shows the potential well of depth Vo and nuclear radius Ro of deuterium 
with binding energy of —2.22 MeV. The next part shows the radius dependence 
of the deuterium radial wave function Xd( r )- The wave-function extends far outside 
the nuclear radius with appreciable amplitude due to the loose binding of deuterium 
ground state. The p-p wave-function Xpp( r ) which comprise the U — initial state 
has small amplitude inside the final nuclear radius. The radial part of the integrand 
entering into the calculation of M space is a convolution of both \d an d Xpp m the 
second and third panels and is given with the hatched shading in the bottom panel. 
It has the major contribution far outside the nuclear radius. 



reaction) as the starting point. The only other reactions with a S(0) greater 
than the above are: d(d, p)t , d(d, n) 3 He , d( 3 He, p) 

4 He, and d( 3 He, -f) 5 Li. However, because of the overwhelmingly large number 
of protons in the stellar thermonuclear reactors, the process involving protons 
on deuterium dominates. The rate of this reaction is so fast compared to its 
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Fig. 6. The Direct Capture reaction d(p, -y) 3 He to form 3 He in its ground state. The 
proton projectile (shown as a plane wave) radiates away a bremsstrahlung photon 
to be captured in a "valence" orbital around the 2 D. 



precursor: p + p — > d + e + v e , that the overall rate of the pp-chain is not 
determined by this reaction. 

One can show that the abundance ratio of deuterium to hydrogen in a 
quasi-equilibrium has an extremely small value, signifying that deuterium 
is destroyed in thermonuclear burning. The time dependence of deuterium 
abundance D is: 

dD 

= t pp - r pd = — < ov > pp -HD < av > pd (28) 

The self regulating system eventually reaches a state of quasi-equilibrium 
and has: 

(D/H) =< av > pp /(2 < av > pd ) = 5.6 x 10" 18 (29) 

at T 6 = 5 and 1.7 x 10~ 18 at T 6 = 40. For the solar system however, this 
ratio is 1.5 x 1CP 4 and the observed (D/H) t, s ratio in the cosmos is ~ 10 -5 . 
The higher cosmic ratio is due to primordial nucleosynthesis in the early phase 
of the universe before the stars formed. (The primordial deuterium abundance 
is a key quantity used to determine the baryon density in the universe). Stars 
only destroy the deuterium in their core due to the above reaction. 
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3.3 3 He burning 

The pp-chain-I is completed (see Fig. 0|) through the burning of 3 He via the 
reaction: 

3 He + 3 Hc^ p + p + 4 He (30) 

with an S-factor: 5(0) = 5500 keV barn and Q-value = 12.86 MeV. In 
addition, the reaction: 

3 He + D -> 4 He +p (31) 

has an S-factor: 5(0) = 6240 keV barn, but since the deuterium concentra- 
tion is very small as argued above, the first reaction dominates the destruction 
of 3 He even though both reactions have comparable 5(0) factors. The cross 
section for the former reaction near the Gamow energy for the sun, has been 
measured in [30] . 

3 He can also be consumed by reactions with He (the latter is pre-existing 
from the gas cloud from which the star formed and is synthesized in the early 
universe and in Pop III objects). These reactions proceed through Direct Cap- 
tures and lead to the ppll and ppIII parts of the chain (happening 15% of 
the time). Note that the reaction 3 He(a, r y) 7 Be together with the subsequent 
reaction: 7 Be(p, 7) 8 £? control the production of high energy neutrinos in the 
sun and are particularly important for the 37 Cl solar neutrino detector con- 
structed by Ray Davis and collaborators. 

3.4 Reactions involving 7 Be 

As shown in Fig.[4l about 15% of the time, 3 He is burned with 4 He radiatively 
to 7 Be. Subsequent reactions involving 7 Be as a first step in alternate ways 
complete the fusion process: 4H — ► A He in the ppll and ppIII chains. 

Electron capture process 

The first step of the ppll chain is the electron capture reaction on 7 Be : 
7 Be + e~ — ► 7 Li + v e (see Fig [7]). This decay goes both to the ground state 
of 7 Li as well as to its first excited state at Ex — 0.478 keV, J 7r = \ ) - 
the percentage of decays to the excited state being 10.4 % in the laboratory. 
The energy released in the reaction with a Q-value of 0.862 keV is carried 
away by escaping mono-energetic neutrinos with energies: E v = 862 and 384 
keV. The measured laboratory mean life of the decay is r = 76. 9d. The 
capture rate in the laboratory can be obtained from Fermi's Golden Rule 
and utilizing the fact that the wave-functions of both the initial nucleus and 
the final one vanish rapidly outside the nuclear domain and the electron wave- 
function in that domain can be approximated as its value at r — and the 
neutrino wave- function by a plane wave normalizes to volume V, so that fljj = 
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& e (0)g/VV J ^* L fi7 Be dT = & e (0)gM n /W, where M n represents the nuclear 
matrix element and the resultant capture rate is: 

Xec = I/tec = (g 2 M*/irc 3 h 4 )E*\%(0)\ 2 (32) 

In the laboratory capture process, any of the various electron shells contribute 
to the capture rate; however the K-shell gives the dominant contribution. At 
temperatures inside the sun, e.g. Tq — 15, nuclei such as 7 Be are largely ion- 
ized. The nuclei however are immersed in a sea of free electrons resulting from 
the ionization process and therefore electron capture from continuum states 
is possible (see e.g., [45], [46]). Since all factors in the capture of continuum 
electrons in the sun are approximately the same as those in the case of atomic 
electron capture, except for the respective electron densities, the 7 Be lifetime 
in a star, t s is related to the terrestrial lifetime Tt by: 

T fr 2|^ t (0)| 2 

rt ~ |JP/r(0)|> (33) 

where |'?/r(0)| 2 is the density of the free electrons n e = p/run at the nucleus, 
p being the stellar density. The factor of 2 in the denominator takes care 
of the two spin states of calculation of the At whereas the corresponding 
Xf r is calculated by averaging over these two orientations. Taking account 
of distortions of the electron wave-functions due to the thermally averaged 
Coulomb interaction with nuclei of charge Z and contribution due to hydrogen 
(of mass fraction Xu) and heavier nuclei, one gets the continuum capture rate 
as: 

= 2|^(0)| 2 r t 

fr ( P /M H )[(l+X H )/2]2nZa(m eC y3kT)i/2 [ > 

with |^ e (0)| 2 ~ (Z/a ) 3 /n. Bahcall et al [44 obtained for the 7 Be nucleus a 
lifetime: 



T fr { 7 Be) = 4.72 x 10 8 



T 



1/2 



p(l+X H ) 

The temperature dependence comes from the nuclear Coulomb field cor- 
rections to the electron wave-function which are thermally averaged. For 
solar condition the above rate |46j gives a continuum capture lifetime of 
Tf r ( 7 Be) = lAOd as compared to the terrestrial mean life of r t = 76.9<i. 
Actually, under stellar conditions, there is a partial contribution from some 
7 Be atoms which are only partially ionized, leaving electrons in the inner K- 
shell. So the contributions of such partially ionized atoms have to be taken 
into account. Under solar conditions the K-shell electrons from partially ion- 
ized atoms give another 21% increase in the total decay rate. Including this, 
gives the solar lifetime of a 7 Be nucleus as: Tq(' Be) = 120d. In addition, the 
solar fusion reactions have to be corrected for plasma electrostatic screening 
enhancement effects. For a recent discussion of the issues see [4"T1 148j . 
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Fig. 7. Electron capture on 7 Be nucleus. The capture proceeds 10.4% of the time 
to the first excited state of 7 Li at 478 keV, followed by a decay to the ground state 
by the emission of a photon. The average energy of the escaping neutrinos (which 
are from the ppll chain) is 814 keV. 

Capture reaction leading to 8 B 

Apart from the electron capture reaction, the Be that is produced is partly 
consumed by proton capture via: 7 Be(p, a) & B reaction. Under solar condi- 
tions, this reaction happens only 0.02% of the time. The proton capture on 
7 Be proceeds at energies away from the 640 keV resonance via the direct cap- 
ture process. Since the product 7 Li nucleus emits an intense 7-ray flux of 478 
keV, this prevents the direct measurement of the direct capture to ground state 
7-ray yield. The process is studied indirectly by either the delayed positron or 
the breakup of the product 8 B nucleus into two alpha particles. This reaction 
has a weighted average S(0) = 0.0238 keVbarn g9]. The 7 Be(p,a) 8 B reac- 
tion cross section measurement has been attempted both by direct capture 
reactions as well as by the Coulomb dissociation of 8 B. For a comparison of 
the 617(0) factors determined by the two methods and a critical review of the 
differences of direct and indirect methods, see [50] . 

The product a B is a radioactive nucleus with a lifetime r = 1.1 s: 



The positron decay of a B(J 7T = 2 + ) goes mainly to the r = 1.6 MeV broad 
excited state in 8 Be at excitation energy E x = 2.94 MeV (J 71 " = 2 + ) due to 
the selection rules (see Fig. [5]). This excited state has very short lifetime and 
quickly decays into two cv-particles. This completes the ppIII part of the pp- 
chain. The average energy of the neutrinos from S B reactions is: E^i^B) = 7.3 
MeV. These neutrinos, having relatively high energy, play an important role 
in several solar neutrino experiments. The neutrino spectrum is not the same 



S B 



8 



Be + e + + v e 



(35) 
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as that obtained from allowed /3-decay approximation and is affected by the 
broad a-unstable 8 Be final state. Winter et al [FT] measured the total energy 
of the a-particles emitted following the /3-decay of 8 B and determined the 
neutrino energy spectrum corrected for recoil order effects in the 8 Be final 
state and constructed from the decay strength function. 




Fig. 8. The decay scheme of S B with positron emission, which goes to the first ex- 
cited state of 8 Be at E x = 2.9 MeV with a width of F = 1.6 MeV. The 8 Be nucleus 
itself fissions into two alpha particles. The neutrinos accompanying the positron 
decay of 8 B are the higher energy solar neutrinos with E v — 7.3 MeV. 



4 The CNO cycle and hot CNO 

The sun gets most of its energy generation through the pp-chain reactions (see 
Fig. [9]). However, as the central temperature (in stars more massive than the 
sun) gets higher, the CNO cycle (see below for reaction sequence) comes to 
dominate over the pp-chain at Tq near 20 (this changeover assumes the solar 
CNO abundance, the transition temperature depends upon CNO abundance 
in the star; in fact, if one is able to isolate the neutrino flux from the Sun's 
weak CN cycle, as in say future extended Borexino or SNO+ experiments 



Thermonuclear fusion and collapse and explosion in massive stars 



25 



with the capability to detect low energy neutrinos, this could in turn directly 
constrain the metallicity of the Sun's corcP^I 38J . The early generation of stars 
(usually referred to as the Population II (Pop II) stars, although there is an 
even earlier generation of Pop III metal poor massive StarJet) generated energy 
primarily through the pp-chain. The Pop II stars are still shining in globular 
clusters, and being of mass lower than that of the sun, are very old. Most other 
stars that we see today are later generation stars formed from the debris of 
heavier stars that contained heavy elements apart from (the most abundant) 
hydrogen. Thus in the second and third generation stars (which are slightly 
heavier than the sun) where higher central temperatures are possible because 
of higher gravity, hydrogen burning can take place through faster chain of 
reactions involving heavy elements C, N, and O which have some reasonable 
abundance (exceeding 1%)) compared to other heavy elements like Li, Be, B 
which are extremely low in abundance. The favored reactions involve heavier 
elements (than those of the pp-chain) which have the smallest Coulomb bar- 
riers but with reasonably high abundance. Even though the Coulomb barriers 
of Li, Be, B are smaller than those of C, N, O (when protons are the lighter 
reactants (projectiles)), they lose out due to their lower abundance. 

In 1937-1938, Bethe and von Weizsacker independently suggested the CN 
part of the cycle, which goes as: 

12 C(p, j) 13 N(e + v e ) 13 C(p, ~/) 14 N(p, j) 15 0(e + v) 15 N(p, a) 12 C (36) 

This has the net result, as before: Ap He + 2e + + 2v e with &Q — 26.73. 
In these reactions, the 12 C and 14 N act merely as catalysts as their nuclei 
are "returned" at the end of the cycle. Therefore the 12 C nuclei act as seeds 
that can be used over and over again, even though the abundance of the 
seed material is minuscule compared to the hydrogen. But note that there is 
a loss of the catalytic material from the CN cycle that takes place through 
the 15 N(p, 7) le O reactions. However, the catalytic material is subsequently 
returned to the CN cycle by the reaction: 16 0(p, -f) 17 F(e + v e ) 17 0(p, a) 14: N. 

In the CN cycle (see Fig fTO]) . the two neutrinos involved in the beta decays 
(of 13 N (t 1/2 = 9.97min) and 15 (t 1/2 = 122.24s)) are of relatively low en- 
ergy and most of the total energy Q = 26.73 MeV from the conversion of four 
protons into helium is deposited in the stellar thermonuclear reactor. The rate 
of the energy production is governed by the slowest thermonuclear reaction in 
the cycle. Here nitrogen isotopes have the highest Coulomb barriers in charged 
particle reactions, because of their Z — 7. Among them 14 N(p, 7) 15 is the 



A measurement CN-cycle neutrino flux (with an expected total flux of about 5 x 
10 s cm -2 s _1 ) would test an assumption of the Standard Solar Model that during 
the early pre-main-sequence Hayashi phase the Sun became homogeneous due to 
convective mixing and that subsequent evolution has not appreciably altered the 
distribution of metals [38] . 

Though subsequently when Carbon has been synthesized by triple-a process, 
these Pop III stars turn on their CN cycle. 
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Fig. 9. Comparison of the temperature dependence of the p-p chain and the CNO 
cycle energy production. The points marked for the solar central temperature T Q = 
Te = 15.7 are shown on both graphs. The CNO cycle generation dominates over 
the pp-chain at temperatures higher than Te = 20, so that for sun like stars, the 
pp-chain dominates. For more massive stars, the CNO cycle dominates as long as 
one of the catalysts: C, N, or O have initial mass concentration at least 1%. Note 
the logarithmic scales of the graph and how both rates drop sharply with decreasing 
temperature, with that of CNO cycle even more drastic due to higher Coulomb 
barriers. 



slowest because this reaction having a final state photon is governed by electro- 
magnetic forces while that involving the other nitrogen isotope: 15 N(p, a) 12 C 
is governed by strong forces and is therefore faster. 

From the CN cycle, there is actually a branching off from 15 N by the reac- 
tion 15 N(p, -f) le O mentioned above. This involves isotopes of oxygen, and is 
called the ON cycle; finally the nitrogen is returned to the CN cycle through 
14 N. Together the CN and the ON cycles, constitutes the CNO bi-cyclc. The 
two cycles differ considerably in their relative cycle-rates: the ON cycle oper- 
ates only once for every 1000 cycles of the main CN cycle. This can be gauged 
from the S(0) factors of the two sets of reactions branching off from 15 N: for 
the 15 N(p,a) 12 C reaction S*(0) = 65 MeV b, whereas for 15 N(p, 7 ) 16 0, it is 
64 keV b, i.e. a factor of 1000 smaller. 
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Fig. 10. The various CNO cycles. The left part is the CN cycle where only C and 
N serve as catalysts for the conversion of four protons into 4 He. Here the slowest 
fusion reaction is (p,7) reaction on 14 N whereas the slower /3-decay has a half-life of 
9.97m. In the CNO bi-cycle (right part), there is leakage from the CN cycle to the 
ON cycle through the branching at 15 N. The flow is returned to the CN cycle (which 
cycles 1000 times for each ON cycle) through 17 0(p, a) 14 N . The right bottom part 
represents additional cycles linking into the CNO cycle through the ir O(p, ~f) 18 F 
reaction [13] , 

4.1 Hot CNO and rp-process 

The above discussion of CNO cycle is relevant for typical temperatures Tq < 
80. These are found in quiescently hydrogen burning stars with solar compo- 
sition which are only slightly more massive than the sun. There are situations 
where the hydrogen burning takes place at temperatures (T ~ 10 8 — 10 9 K) 
which are in excess of those found in the interiors of the ordinary "main 
sequence" stars. Examples of these are: hydrogen burning at the accreting 
surface of a neutron star or in the explosive burning on the surface of a white 
dwarf, i.e. novae, or the outer layers of a supernova shock heated material in 
the stellar mantle. These hot CNO cycles operate under such conditions on 
a rapid enough time scale (few seconds) so that even "normally" /3-unstable 
nuclei like 13 N will live long enough to be burned by thermonuclear charged 
particle reactions, before they are able to /3-decay [55], [53]. So, unlike the 
normal CNO the amount of hydrogen to helium conversion in hot CNO is 



28 A. Ray 



limited by the /3-decay lifetimes of the proton-rich nuclei like: 14 and lb O 
rather than the proton capture rate of 14 N, Wallace and Woosley [51] showed 
that for temperatures, T > 5 x lO 8 !^ , nucleosynthesised material can leak out 
of the cycles. This leads to a diversion from lighter to heavier nuclei and is 
known as the rapid proton capture or rp-process. The flow between the hot 
CNO cycle and the rp capture process in X-ray bursts from the atmosphere 
of a neutron star is regulated by the 15 0(a,j) Ne reaction. Another alpha 
capture reaction ls Ne(a, 7) 21 iVa, continuously processes the available A He 
nuclei flowing towards heavier elements. For a recent discussion on the hot 
CNO and the rp process on accreting neutron stars, see [55] , 

The nucleosynthesis path of rp-process of rapid proton addition is anal- 
ogous to the r-process of neutron addition (for neutron capture processes in 
the early galaxy see [56]). The hot hydrogen bath converts CNO nuclei into 
isotopes near the region of proton unbound nuclei (the proton drip line). 
For each neutron number, a maximum mass number A is reached where the 
proton capture must wait until /3 + -decay takes place before the buildup of 
heavier nuclei (for an increased neutron number) can take place. Unlike the 
r-process the rate of the rp-process is increasingly hindered due to the in- 
creasing Coulomb barrier of heavier and higher-Z nuclei to proton projectiles. 
Thus the rp-process does not extend all the way to the proton drip line but 
runs close to the beta-stability valley and runs through where the /3 + -decay 
rate compares favorably with the proton captures. 

5 Helium burning and the triple-e* reaction 

After hydrogen burning in the core of the star has exhausted its fuel, the 
helium core contracts slowly. Its density and temperature goes up as gravita- 
tional energy released is converted to internal kinetic energy. The contraction 
also heats hydrogen at the edge of the helium core, igniting the hydrogen to 
burn in a shell. At a still later stage in the star's evolution, the core has con- 
tracted enough to reach central temperature density conditions: Tq = 100—200 
and p c = 10 2 — 10 5 gm cm -3 when the stellar core settles down to burn A He 
in a stable manner. The product of helium burning is 12 C. Since in nature, 
the A = 5 and A = 8 nuclei arc not stable, the question arises as to how 
helium burning bridges this gap. A direct interaction of three a particles to 
produce a 12 C nucleus would seem at first sight, to be too improbable (as 
was mentioned, for example, in Bethe's 1939 paper [57], which was primarily 
on the CN cycle). However, Opik [58] and Salpeter [43], [59] independently 
proposed a two step process where in the first step, two a particles interact to 
produce a s Be nucleus in its ground state (which is unstable to a-breakup), 
followed by the unstable nucleus interacting with another a-particle process 
to produce a 12 C nucleus. 

Thus the triple alpha reaction begins with the formation of 8 Be that has 
a lifetime of only 1 x 10~ 16 s (this is found from the width r = 6.8 eV of the 
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ground state and is the cause of the A — 8 mass gap). This is however long 
compared to the transit time 1 x 10~ 19 s of two a-particles to scatter past 
each other non-resonantly with kinetic energies comparable to the Q-value of 
the reaction namely, Q = —92.1 keV. So it is possible to have an equilibrium 
build-up of a small quantity of s Be in equilibrium with its decay or reaction 
products: a + a Be. The equilibrium concentration of the s Be nucleus 
can be calculated through the Saha equation 

N 12 = (2.7+1) ( _Ek } (3?) 

12 2 y nkT> (2Ji + 1)(2J 2 + 1) kT y V ' 

at the relevant temperature Tq = 11 and p = 10 5 gm cm -3 to be: 

Salpeter suggested that this small quantity of 8 Be serves as the seed for 
the second stage of the triple a-capture into the 12 C nucleus. It was however 
shown by Hoyle [HI] that the amount of 12 C produced for the conditions inside 
a star at the tip of the red-giant branch is insufficient to explain the observed 
abundance, unless the reaction proceeds through a resonance process [60] . 
The presence of such a resonance greatly speeds up the rate of the triple-a 
process which then proceeds through an s-wave (I = 0) resonance in 12 C near 
the threshold of 8 Be + a reaction. Since 8 Be and 4 He both have J 71 ' = + , 
an s-wave resonance would imply that the resonant state in question has 
to be + in the 12 C nucleus. Hoyle suggested the excitation energy to be: 
Ex ~ 7.68 MeV in the 12 C nucleus and this state was experimentally found 
by W.A. Fowler's group ([62]) with spin-parity: J 71 — + . This state has 
a total width ( 13J ) r = 8.9 ± 1.08 eV, most of which lies in r a , due to the 
major propensity of the 12 C nucleus to break-up through a-decay. (The decay 
of the excited state of 12 C by 7-rays cannot go directly to the ground state, 
since the resonance state as well as the ground state of the 12 C nucleus have 
both J 71 " = + and + — * + decays are forbidden. This partial width due to 
gamma-decay is several thousand times smaller than that due to a-decay). 
So, r — r a + r ra d ~ r a and r rac i = r i + r e + e - = 3.67 ± 0.50 meV. Again 
the radiative width r ra( i is dominated by the width due to photon width 
deexcitation: = 3.58 ± 0.46 meV. (Note the scales of mi'&lectron Volts). 
The reaction scheme for the first and the second parts of the triple-alpha 
reaction is given in Fig. 1111 The locations of the Gamow energy regions near 
the above resonance state (for several stellar temperatures) are shown only 
schematically. 

The reaction rate for the 12 C formation can be calculated by using the 
properties of the resonant state and the thermally averaged cross-section: 

r3a = Ns Be N a < av >s Be+a (39) 
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Fig. 11. Triple alpha process of C synthesis. In the first step a small amount of 
8 Be nuclei builds up in equilibrium with its decay products (forward and backward 
reactions involve alpha particles). The second step involves a capture of another 
alpha particle by the unstable 8 Be nucleus which proceeds via an s-wave resonance 
state in the product nucleus 12 C located close to the Gamow energy window for 
temperatures indicated schematically by the three-way arrows on the right. 



Here Ns Be and N a are the number densities of interacting s Be and 4 He 
nuclei and the angular brackets denote thermal averaging over a Maxwell 
Boltzmann distribution tp(E). This averaging leads to: 



r 3a = Ns Be N a f ^(E)v(E)a(E)dE (40) 
Jo 



with 
and 

a(E) = 



2 E c\F, 



V5rfcT v ' ; (kTE)V2 

a N 2 2j + i r t r 2 



^2ir> (2Ji + 1)(2J 2 + 1) (E - E R ) 2 + (r/2) 2 

is the Breit-Wigner resonant reaction cross section with the resonant energy 
centroid at E = Er. The total width r is a sum of all decay channel widths 
such as i~i = r a and _T 2 = -T 7 - If the width r is only a few eVs then the 
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functions ip(E) and v(E) can be pulled out of the integral. Then, the reaction 
rate will contain an integral like: J Q aB\v(E)dE — 2Tr(\/2irfr) 2 ujrir2/r, 
where w = (2 J + l)/[(2Ji + 1)(2J 2 + 1)] and the functions pulled out of the 
integral need to be evaluated at E = Er. Since most of the time the excited 
state of the 12 C* breaks-up into a-particles, we have A = r a dominating 
over _T 7 and {Txr^/r) ~ i~2. This limit usually holds for resonances of energy 
sufficiently high so that the incident particle width ) dominates the natural 
width of the state (I^)- In that case, we can use the number density of the 
s Be nuclei in equilibrium with the a-particle nuclei bath as described by Saha 
equilibrium condition: 

NCBe) = Ar> / __^l_exp(-E r /kT) (41) 

where f is the screening factor. It is possible to get the overall triple- 
alpha reaction rate by calculating the equilibrium concentration of the excited 
(resonant) state of 12 C reached by the 8 -Be + a -^ 12 C* reaction and then 
multiplying that concentration by the gamma-decay rate r^/Ti which leads 
to the final product of 12 C. So, the reaction rate for the final step of the 
triple-alpha reaction turns out to be: 



r 3a =N SBe N a K 2 (^=) w/r 2 exp(-E r /kT) (42) 



3/2 

\We) 

where fi is the reduced mass of the reactants s Be and a particle. This 
further reduces by the above argument to: 

rs „„ c .f^(^l) 3 /Vl« P (- Q / k T, (43, 

The Q-value of the reaction is the sum of Eft( 8 Be + a) = 287 keV and 
E R (a + a) = \Q\ = 92 keV and turns out to be: Q 3a = (Mi2 C * - 3M Q )c 2 = 
379.38 ± 0.20keV [53]. Numerically, the energy generation rate for the triple- 
alpha reaction is: 

e 3a = = 3.9 x lO 11 ^/ exp(~42.94/T 8 ) erg gnr 1 s" 1 (44) 

The triple alpha reaction has a very strong temperature dependence: near 
a value of temperature To, one can show that the energy generation rate is: 

e(T) = e(T )(^r (45) 

where, n = 42.9/Ts — 3. Thus at a sufficiently high temperature and den- 
sity, the helium gas is very highly explosive, so that a small temperature rise 
gives rise to greatly accelerated reaction rate and energy liberation. When 
helium thermonuclear burning is ignited in the stellar core under degenerate 
conditions, an unstable and sometimes an explosive condition develops. 
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6 Survival of 12 C in red giant stars and 12 C(ct,~f) 16 
reaction 

The product of the triple-alpha reactions 12 C, is burned into le O by a-capture 
reactions: 

12 C* + a-> 16 + 7 (46) 

If this reaction proceeds too efficiently in helium burning Red giant stars, 
then all the carbon will be burned up to oxygen. Carbon is however the most 
abundant element in the universe after hydrogen, helium and oxygen, and the 
cosmic C/0 ratio is about 0.6. In fact, the O and C burning reactions and the 
conversion of He into C and O take place in similar stellar core temperature 
and density conditions. Major ashes of He burning in Red Giant stars are 
C and O. Red Giants are the source of the galactic supply of 12 C and le O. 
Fortuitous circumstances of the energy level structures of these alpha-particle 
nuclei are in fact important for the observed abundance of oxygen and carbon. 

For example, if as in the case of the 3a reaction, there was a resonance in 
the 12 C(a,7) 16 reaction near the Gamow window for He burning conditions 
(i.e. T 9 ~ 0.1 — 0.2), then the conversion of 12 C — ► le O would have proceeded 
at a very rapid rate. However, the energy level diagram of le O shows that for 
temperatures upto about Tg ~ 2, there is no level available in 16 O to foster 
a resonant reaction behavior (Fig. I12p . But since this nucleus is found in na- 
ture, its production must go through either: 1) a non-resonant direct capture 
reaction or 2) non-resonant captures into the tails of nearby resonances (i.e. 
sub-threshold reactions). In Fig.021 also shown on the left of the 16 O energy 
levels, is the threshold for the 12 C + 4 He reaction, drawn at the appropriate 
level with respect to the ground state of the le O nucleus. The Gamow energy 
regions drawn on the extreme right for temperatures Tg = 0.1 and above, in- 
dicates that for the expected central temperatures, the effective stellar (center 
of mass) energy region is near Eq = 0.3 McV. This energy region is reached 
by the low energy tail of a broad resonance centered at Ecu = 2.42 MeV 
above the threshold (the J 71 = 1~ state at 9.58 MeV above the ground state 
of 16 O) with a (relatively large) resonance width of 400 keV. On the other 
hand, there are two sub-threshold resonances in le O (at Ex = 7.12 MeV and 
Ex = 6.92 MeV), i.e. -45 keV and -245 keV below the a-particle threshold 
that have J T = l - and J w = 2 + , that contribute to stellar burning rate 
by their high energy tails. However, electric dipole (El) 7-decay of the 7.12 
MeV state is inhibited by isospin selection rules. Had this not been the case, 
the 12 C(a,j) 16 reaction would have proceeded fast and 12 C would have 
been consumed during helium burning itself. The two sub-threshold states 
at —45 keV and —245 keV give contributions to the astrophysical S-factor of: 
Sx- (E ) = 0.1 MeV barn and S 2 + (E ) = 0.2 MeV barn respectively at the rel- 
evant stellar energy E = 0.3 MeV. The state at E C m = 2.42 MeV (J* = 1~ 
state at 9.58 MeV) gives a contribution: Sx- (E Q ) = 1.5 x 10~ 3 MeV barn. The 



Thermonuclear fusion and collapse and explosion in massive stars 



33 



total S-factor at E = 0.3 MeV is therefore close to 0.3 MeV barn. These then 
provide low enough S or cross-section to not burn away the 12 C entirely to 
16 0, so that C/O ~ 0.1 at the least. 
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Fig. 12. Energy levels of O nucleus near and above the alpha-particle threshold of 
capture on 12 C. Shown on the right are effective stellar energy regions corresponding 
to the temperatures given near the three-way arrows. The reaction rate is influenced 
mainly by the high energy tails of two sub-threshold resonances in 16 O at Er = 
—45 keV and Er — — 245 keV, plus the low energy tail of another high-lying broad 
resonance at 9580 keV. 



Additionally, le O nuclei are not burnt away by further a-capture in the 
reaction: 

le O+ 4 He^ 20 iVe + 7 (47) 

A look at the level schemes of 20 Ne (see Fig. IT3"]) shows the existence of 
a Ex = 4.97 MeV state (J n = 2~) in the Gamow window. However, this 
state cannot form in the resonance reaction due to considerations of parity 
conservation (unnatural parity of the resonant statef^l The lower 4.25 MeV 
state (J 77 = 4 + ) in 20 Ne also cannot act as a sub-threshold resonance as it 

17 Whether or not a resonant state can be formed or accessed via a given reaction 
channel depends upon the angular momentum and parity conservation laws. The 
spins of the particles in the entrance channel, and relative angular momen- 
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lies too far below threshold and is formed in the g-wave state. Therefore only 
direct capture reactions seem to be operative, which for (a, 7) reactions lead 
to cross-sections in the range of nanobarns or below. Thus the destruction of 
le O via: 16 0(a, 7) 20 Ne reaction proceeds at a very slow rate during the stage 
of helium burning in Red Giant stars, for which the major ashes are carbon 
and oxygen and these elements have their galactic origin in the Red Giants. 

To summarize, the synthesis of two important elements for the evolution 
of life as we know on the earth have depended upon fortuitous circumstances 
of nuclear properties and selection rules for nuclear reactions. These are: 1) 
the mass of the unstable lowest (ground) state of 8 Be being close to the 
combined mass of two a-particles; 2) there is a resonance in 12 C at 7.65 MeV 
which enhances the alpha addition reaction (the second step); and 3) parity 
conservation has protected le O from being destroyed in the le O(a, j) 20 Ne 
reactions by making the 4.97 MeV excited state in 20 Ne of unnatural parity. 

The experimental determination of the reaction rate 12 C(a, 7) le O has been 
an important goal in nuclear astrophysics for several decades. Its cross sec- 
tion at the position of the Gamow window for a typical stellar temperature of 
2.5 x 10 8 K is comparable to that of weak interaction cross-sections. At those 
energies, this reaction is practically a non-resonant reaction and its cross- 
section is determined by the tails of interfering resonance and sub-threshold 
states [64] . The low cross section and the complexity of low energy contribu- 
tions to the reaction rate makes a reliable prediction difficult [65j [66] • 

7 Advanced stages of thermonuclear burning 

As the helium burning progresses, the stellar core is increasingly made up of 
C and O. At the end of helium burning, all hydrogen and helium is converted 
into a niixtur^f] of C and O, and since H, He are most abundant elements 
in the original gas from which the star formed, the amount of C and O are 
far more in the core than the traces of heavy elements in the gas cloud. 
Between these two products, the Coulomb barrier for further thermonuclear 
reaction involving the products is lower for C nuclei. At first the C+O rich 



turn I adds upto the angular momentum of the resonant state J = ji+js+l. There- 
fore, for spin-less particles like the closed shell nuclei 4 He, 1(i O (ji = 0, j'2 = 0), 
we have J = Z. In the entrance channel of the reacting particles, the parity would 
be: (— 1) 7r(ji)7r(j2) = ( — l) i=0 (1) (1) . If the parity of the resonance state were 
the same as that of the entrance channel, then the resultant state would have 
been a "natural parity" state. However, since the 4.97 MeV state in 20 Ne has an 
assignment: J" = 2~, this is an "unnatural parity" state. 

Note however the caveat: if the amount of 12 C is little (either due to a long stellar 
lifetime of He burning or due to a larger rate of the 12 C + a O + 7 reaction 
whose estimate outlined in the earlier section is somewhat uncertain), then the 
star may directly go from He-burning stage to the O-burning or Ne-burning stage 
skipping C-burning altogether f|67]h 
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Fig. 13. Energy levels of nuclei participating in thermonuclear reactions during the 
helium burning stage in red giant stars (after [13]). Survival of both 12 C and le O 
in red giants, from which terrestrial abundances result, depends upon the fortuitous 
circumstances of nuclear level structures and other properties of these nuclei. 



core is surrounded by He burning shells and a helium rich layer, which in turn 
may be surrounded by hydrogen burning shell and the unignited hydrogen 
rich envelope. When the helium burning ceases to provide sufficient power, 
the star begins to contract again under its own gravity and as implied by 
the Virial theorem the temperature of the helium exhausted core rises. The 
contraction continues until either the next nuclear fuel begins to burn at rapid 
enough rate or until electron degeneracy pressure halts the infall. 

7.1 Carbon burning 

Stars somewhat more massive than about 3 M Q contract until the temperature 
is large enough for carbon to interact with itself (stars less massive on the main 
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sequence may settle as degenerate helium white dwarfs). For stars which are 
more massive than M > 8— 10 M Q (mass on the main sequence, - not the mass 
of the C+O core), the contracting C+O core remains non-degenerate until C 
starts burning at T ~ 5 x 10 8 K and p = 3 x 10 6 gcm~ 3 . Thereafter sufficient 
power is generated and the contraction stops and quiescent (hydrostatic, not 
explosive) C-burning proceeds (see Fig. [T4|) . 




Fig. 14. Tracks in the core temperature, density plane of stars of various masses 
(at the start of hydrogen burning i.e. main sequence masses). Note that a star of 
mass M ~ 15Mq ignites all its fuels in non-degenerate conditions, whereas a star of 
mass M ~ 4Mq ignites carbon under strongly degenerate conditions. (After [68]). 

The combined mass of two reacting 12 C nuclei falls at an excitation energy 
of 14 MeV in the compound nucleus of 24 Mg. At this energy there are many 
compound nuclear states, and the most effective range of stellar energies (the 
Gamow window) at the relevant temperature is about 1 MeV; hence a number 
of resonant states can contribute to the decay of the compound nucleus, and 
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even the large angular momentum resonances may be important because the 
penetration factors in massive nuclei are not affected by centrifugal barriers. 
The carbon on carbon burning can proceed through multiple, energetically 
allowed reaction channels, listed below: 

12 C+ 12 C^ 20 Ne+ 4 He (Q = 4.62 MeV) (48) 
-> 23 Na+p (Q = 2.24 MeV) (49) 

-> 23 Mg + n (Q = -2.62 MeV) (50) 

At the temperatures where carbon burning starts, the neutron liberating 
reactions requires too much particle kinetic energy to be effective. In addition, 
based on laboratory measurements at higher energies compared to the stellar 
energies, the electromagnetic decay channel ( 24 Mg + 7) and the three parti- 
cle channel ( le O + 2cv) have lower probability compared to the two particle 
channels: 23 Na + p and 20 Ne + a. The latter two channels have nearly equal 
probabilities (see [TJ; at the lowest center of mass energies for which cross- 
sections are measured in the laboratory for the proton and a channels, (i.e. 
about 2.45 MeV [69]), the branching ratios were b p ~ 0.6 and b a ~ 0.4), and 
therefore the direct products of carbon burning are likely to be 23 Na, 20 Ne, 
protons and alpha particles. The rate for this reaction per pair of 12 C nuclei 
is ([70]): 

logA 12 , 12 = logf 12 , 12 + 4.3 - 36.55(1 + 0TT 9 )V3 _ 2 ^ 

the factor /i2,i2 is a screening factor. Now, at the temperatures of 12 C 
burning, the liberated protons and alpha particles can be quickly consumed 
through the reaction chain: 12 C(p, j) 13 N(e + v e ) 13 C(a, n) 16 0. Thus, the net 
effect is that the free proton is converted into a free neutron (which may 
be further captured) and the a-particle is consumed with 12 C into 16 0. The 
a-particles are also captured by other alpha-particle nuclei, resulting in, at 
the end of carbon burning in nuclei like: 16 0, 20 Ne, 2i Mg and 28 Si. These 
secondary reactions augment the energy released by the initial carbon reaction 
and Reeves (1959) estimated that each pair of 12 C nuclei release about 13 
MeV of energy. Towards the end of carbon burning phase there are also other 
reactions such as: 12 C + 16 O and 12 C + 20 Ne which take place. But these are 
less rapid and are not expected to play major roles compared to the 12 C+ 12 C 
reactions, due to their increased Coulomb barriers. A recent discussion of the 
heavy ion reactions involving C and O is contained in [T] section 3.6. 

During the carbon-burning and subsequent stages, the dominant energy 
loss from the star is due to neutrinos streaming out directly from the stellar 
thermonuclear furnace, rather than by photons from the surface. The neutrino 
luminosity is a sensitive function of core temperature and quickly outshines the 
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surface photon luminosity of the star at carbon burning stage. The (thermal) 
evolutionary time scale of the star, due to the neutrino emission becomes very 
short and the core evolves rapidly, - so rapidly (compared to the "cooling" 
time scale Kelvin- Hclmholtz time: Tkb ~ GM 2 / RL p h) that the conditions in 
the core are "uncommunicated" to the surface, since this happens by photon 
diffusion. The surface conditions (e.g. the temperature) then does not change 
much as the core goes beyond the carbon burning stage, and it may not be 
possible just by looking at a star's surface conditions whether the core is 
close to a supernova stage or has many thousands of years of hydrostatic 
thermonuclear burning to go. 



7.2 Neon burning 

The result of carbon burning is mainly neon, sodium and magnesium, but 
aluminum and silicon are also produced in small quantities by the capture of 
a, p and n released during carbon burning. When carbon fuel is exhausted, 
again the core contracts and its temperature T c goes up. At approximately 
Tg ~ 1, energetic photons from the high energy tail of the Planck distribution 
function can begin to disintegrate the 20 Ne ash (see Fig. [T3|) so that one has 
the reaction: 20 Ne + 7 O + 4 He. 

Nucleons in a nucleus are bound with typical binding energy of several to 
8 MeV. An energetic 7-ray photon is required to photo-eject a nucleon. Two 
nucleon ejection requires more energy. Alpha particles are however released at 
approximately the same energy as a nucleon due to the low separation energy 
of an alpha particle in the nucleus. For example, the alpha separation energy 
in 20 Ne is 4.73 MeV. Thus, the major photo-nuclear reactions are: (7, n), (7,73) 
and (7, a) processes. For a photo-disintegration reaction to proceed through 
an excited state Ex in the mother, the decay rate is:- 



A(7,a) 



/ 2Jr + i r 7 

6XP ' ~ ~kT' 2Jn + 1 f 



x ^ (52) 



In the above equation, the first factor in square brackets on the RHS is the 
probability of finding the nucleus in the excited state Ex and spin Jr (with 
Jo being the ground state spin), while the second factor r a /Ti is the decay rate 
of the excited state with an alpha particle emission. Now since Ex = En. + Q, 
we have: 

A(7,°) = ^JoTi? {2JR + ^^jfi^-ER/kT) (53) 

At Tg > 1, the photo-disintegration is dominated by the 5.63 MeV level 
in 20 Ne (see Fig. [13]). At approximately Tg ~ 1.5, the photo-dissociation rate 
becomes greater than the rate for alpha capture on 16 to produce 20 A^e (i.e. 
the reverse reaction), thus leading effectively to the net dissociation of 20 Ne. 
The released 4 He reacts with the unspent 20 JVe and leads to: 4 He + 20 Ne — > 
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2i Mg + 7. Thus the net result of the photo-dissociation of two 20 Ne nuclei 
is: 2 x 20 Ne -> 16 + 24 Mg with a net Q-value of 4.58 MeV. The brief neon 
burning phase concludes at Tg close to ~ 1.5. 

7.3 Oxygen burning 

At the end of the neon burning the core is left with a mixture of alpha particle 
nuclei: le O and 2i Mg. After this another core contraction phase ensues and 
the core heats up, until at Tg ~ 2, 16 (9 begins to react with itself: 



The first reaction takes place approximately 45% of the time with a De- 
value of 9.593 MeV. In addition to Si and S, the oxygen burning phase also 
produces, Ar, Ca and trace amounts of CI, K, etc upto Sc. Then at T 9 ~ 3, 
the produced 2S Si begins to burn in what is known as the Si burning phase. 

7.4 Silicon burning 

As we have seen, most of the stages of stellar burning involve thermonuclear 
fusion of nuclei to produce higher Z and A nuclei. The first exception to 
this is neon burning where the photon field is sufficiently energetic to photo- 
dissociate neon, before the temperature rises sufficiently to allow fusion reac- 
tions among oxygen nuclei to overcome their Coulomb repulsion. Processing 
in the neon burning phase takes place with the addition of helium nuclei to 
the undissociated neon rather than overcoming the Coulomb barrier of two 
neon nuclei. This trend continues in the silicon burning phase. In general, a 
photo-disintegration channel becomes important when the temperature rises 
to the point that the Q-value, i.e. the energy difference between the fuel and 
the products is smaller than approximately 30fcsT f |71j). 

With typical Q- values for reactions among stable nuclei above silicon be- 
ing 8-12 MeV, photo-disintegration of the nuclear products of neon and oxy- 
gen burning begins to play an important role once the temperature exceeds: 
Tg > 3. Then nuclei with smaller binding energies are destroyed by photo- 
dissociation in favor of their more more tightly bound neighbors, and many 
nuclear reactions involving a-particles, protons and neutrons interacting with 
all the nuclei in the mass range A = 28 — 65 take place. In contrast to the 
previous burning stages where only a few nuclei underwent thermonuclear 
reactions upon themselves, here the nuclear reactions are primarily of a re- 
arrangement type, in which a particle is photo-ejected from one nucleus and 
captured by another and a given fuel nucleus is linked to a product nucleus 
by a multitude of reaction chains and cycles and it is necessary to keep track 
of many more nuclei (and many reaction processes involving these) than for 



16 +16 O _2S Si + 4 H( , 



(54) 




(55) 
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previous burning stages. More and more stable nuclei form in a nuclear reac- 
tion network as the rearrangement proceeds. Since there exists a maximum in 
the binding energy per nucleon at the 56 Fe nucleus, the rearrangements lead 
to nuclei in the vicinity of this nucleus (i.e. iron-group nuclei). 

In the mass range A = 28 — 65, the levels in the compound nuclei that form 
during silicon burning are so dense that they overlap. Moreover, at the high 
temperatures that are involved (Tg = 3 — 5), the net reaction flux may be small 
compared to the large forward and backward reactions involving a particular 
nucleus and a quasi-equilibrium may ensue between groups of nuclei which 
are connected between separate groups by a few, slow, rate-limiting reactions 
("bottlenecks"). However, as the available nuclear fuel(s) are consumed and 
thermal energy is removed due to escaping neutrinos, various nuclear reactions 
may no longer occur substantially rapidly ("freeze-out"). Thielemann and Ar- 
nett [72] found that for cores of massive stars in hydrostatic cases, the bot- 
tlenecks between quasi-equilibrium (QSE) groups coincided with Z=21 nuclei 
whereas for lower mass stars, lower temperatures and Y e and higher density 
this bridge involved neutron rich isotopes of Ca capturing protons. Hix and 
Thielemann [71] discussed and contrasted these results with those of earlier 
workers; in general the reaction flow across the boundary of the QSE groups 
arc influenced by the neutralization of the material, i.e. the overall Y e . It is 
in this context that weak interaction processes such as electron capture and 
beta decay of nuclei are important, by influencing the Y e and thereby the re- 
action flow. These ultimately affect both the stellar core density and entropy 
structures, and it is important to track and include the changing Y e of the 
core material not only in the silicon burning phase, but even from earlier oxy- 
gen burning phases. The calculation of stellar weak processes on nuclei has 
spawned extensive literature (see [73], [74], [75] etc., and for a review [76]). 

Iron, nickel and neighboring nuclei which are the products of the hydro- 
static Si burning in the core are mostly trapped in the collapsing core that 
ends up as the compact remnant and little of this may reach the interstellar 
medium. However, the blast wave shock launched after the core bounce im- 
pacts through the onion-like outer shells that the star's earlier evolution has 
left behind. The nearest part of the Si shell that is still unburned is heated to 
sufficient temperatures to form iron peak nuclei. Regions which are not as close 
can undergo incomplete Si burning and be left with substantial amounts of Si, 
S, Ar, Ca and Ti. Three separate outcomes may result depending upon the ini- 
tial density and the peak temperature: 1) incomplete Si burning; 2) a "normal 
freezeout" and 3) an "a- rich freezeout" . In the first case, with initial tem- 
peratures of about 5 x 10 9 K and peak density of p = 10 9 g cm -3 , significant 
amounts of Si and other intermediate-mass elements remain after the charged 
particle reactions freezeout in the expanding ejecta. In the normal freezeout, 
an initial condition of 7 x 10 9 K and peak density of p = 10 9 g cm -3 leads to 
complete Si burning. The "a-rich freezeout" however takes place at lower peak 
densities p = 10 7 g cm -3 though similar peak temperatures: 7 x 10 9 K. This 
ends up with an abundance of 4 He nuclei which produces a significant flow 
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through the triple alpha reaction into intermediate mass nuclei, in addition 
to the iron group nuclei products. The feedback between the rate of nuclear 
recombination and that of temperature evolution in the expansion critically 
controls the "cv-richness" of matter and production of important radioactive 
nuclei e.g. 44 Ti and is one of the challenges of computational simulation of 
silicon burning. For a recent discussion of the computational aspects of the 
nuclear evolution during silicon burning, see |77] . 

In summary, a few key points concerning the thermonuclear burning of 
28 Si are as follows :- 

• Direct thermonuclear fusion of two 28 Si nuclei does not take place be- 
cause their Coulomb barrier is too high. Instead thermonuclear fusion takes 
place by successive additions of a-particles, neutrons and protons. 

• Although this is actually a large network of nuclear reactions it is 
called "silicon burning" because 28 Si offers the largest resistance to photo- 
dissociation because of its highest binding energy among intermediate mass 
nuclei. 

• The source of the a-particles which are captured by 2 Si and higher 
nuclei is 28 Si itself. Silicon, sulphur etc. partially melt-down into a-particles, 
neutrons and protons by photo-dissociation. These then participate in reaction 
networks involving quasi-equilibrium clusters linked by "bottleneck" links. 

• Although beta decay and electron captures on stellar core nuclei do not 
produce energy in major ways they nevertheless play a crucial role in shifting 
the pathways of nuclear and thermodynamic evolution in the core conditions. 
These ultimately determine the mass of the core and its entropy structure 
which finally collapses in a supernova explosion. 

8 Core collapse SNe: electron capture and neutrinos 

At the end of nuclear burning in the core of a massive star, inert Iron group 
nuclei are left in the innermost region, - inert because nuclear burning cannot 
extract any further energy from the most tightly bound nuclei. At the last 
stages of nuclear burning, neutrino cooling which is far more efficient than 
the radiation of photons from the surface, would leave the core compact and 
degenerate [78j (where electron Fermi energy is much larger than the thermal 
energy kT) due to a number of processes such as e + /e~ pair annihilation 
process, photoneutrino process and plasmon decay into neutrinos. These neu- 
trinos escape from the star freely at this stage since their interaction cross 
section is so small, much smaller than that of the photons. In the picture 
of the late stages of a spherical star that B 2 FH [79] proposed one had an 
onion-skin model: nested shells of progressively heavier elements, with the 
densest iron core at the center. Successively high-Z elements are ignited in 
increasingly central parts of the core in decreasingly lower specific entropy 
environments, and are contained within entropy barriers left behind by pre- 
vious generations of nuclear burning. The degenerate core becomes unstable 
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to gravitational collapse at roughly the same Chandrasekhar mass (modulo 
relatively small variations in the electron fraction or entropy profiles in the 
core) practically independent of the total mass of the star, which for example 
could range from ~ 8M to ~ 60M Q [HZ] . This "core convergence" establishes 
the important connection between the Chandrasekhar mass and the masses 
of neutron stars that are formed from core-collapse supernovae (though the 
details of the boundaries of the main sequence masses that lead to neutron 
stars, as well as the composition of the cores, e.g. O-Mg-Ne cores vs iron 
cores etc., have evolved as further research was undertaken). The supernova 
explosion itself is triggered by the collapse of the core [80], followed by the 
sudden release of gravitational energy. When the iron core grows in size, its 
temperature increases, and when T > 7x lO 9 ^, iron nuclei photo-dissociate in 
cndothermic reactions (consuming energy, even though the free energy F = U 
-TS decreases due to increased entropy) into alphas and neutrons. This leads 
to the effective adiabatic index decrease below 4/3 and the iron core becomes 
unstable to collapse. Collapse initiation through photo-dissociation happens 
in relatively massive cores of more massive stars whereas in less massive stars 
the collapse may initiate due to reduction of electrons (which provide the bulk 
of the supporting pressure due to their high Fermi momenta) as they are cap- 
tured by nuclei in "inverse beta decay" in regions where the electron Fermi 
energy exceed the capture thresholds. 

The core collapses on a dynamical scale until the infall is suddenly halted 
when the central density overshoots that of the nuclear matter, at 4 — 5 x 
10 14 gm cm" 3 . In between, the iron core collapsed onto itself nearly freely at 
about a quarter of the speed of light. Initially it had the size of the Earth, 
but towards the end it is a hot, dense, neutron rich sphere about 30 km in 
radius (see [81] for a review). A static accretion shock wave forms initially at 
the edge of the quasi-free falling core, which soon starts propagating outward 
through the outer core and mantle as more kinetic energy is brought into it 
by infalling matter. The shock wave soon stalls (at least for most of the range 
of main sequence masses beyond about 10 Mq) which may however be helped 
after about 500 milliseconds ("after a pause that refreshes" - according to 
Hans Bethe) |82] by neutrinos which are freely streaming out from the inner 
core, but manage to deposit enough energy (albeit a small fraction of their 
total energy) "reviving" the shock. During this post-bounce phase the proto- 
neutron star radiates away most of its energy (3 x 10 53 erg) in the form of 
neutrinos and antineutrinos of various flavors and this accounts for the larger 
intrinsic energy per unit mass of the central engines of core collapse SNe. 

8.1 Electron capture on nuclei and protons: a core thermometer 

During the gravitational collapse, the entropy of the core stays low, which 
permit the nuclei of various elements present in the core to (largely) sur- 
vive thermal disintegration and coexist with a small fraction of "dripped" 
nucleons (|83J hereafter BBAL). Around the density 5 x 10 n gm cm" 3 , neu- 
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trinos produced through electron capture on these nuclei and free protons are 
trapped in the core. Much of the information pertaining to the conditions in 
which the neutrinos are originally produced, such as the nuclear and ther- 
modynamic properties of the core of the supernova are altered because these 
neutrinos undergo inelastic scattering with the overlying stellar matter in the 
post neutrino trapping phase. Neutrinos which are emitted through electron 
captures on the nuclei present in the pre-supernova and collapsing core be- 
fore it reaches neutrino trapping density (~ 10 12 gm cm -3 ) [84], [85], however, 
stream freely through the stellar matter without any further interactions. 
These pre-trapping neutrinos carry with them information on both physical 
conditions within the core, as well as it's nuclear configuration e.g. the ratio of 
the number density of free protons to that of heavy nuclei. The last quantity 
can depend on the nuclear equation of state relevant to collapse. Since neu- 
trinos act as probes of the dynamic, thermodynamic and nuclear properties 
of the pre-supernova and collapsing core, their detection and measurement of 
energy spectrum can have significant implications. The time evolution of the 
detected spectrum could also reveal the dynamical time scale - a clue to the 
average density and mass of the stellar core which may have implications for 
neutron star vs black hole formation^!. The reduction of lepton fraction dur- 
ing stellar collapse has implications for shock formation stages and the overall 
dynamics - even in the delayed explosion stage, since it determines, through 
the original energy of the bounce shock, and the entropy profile in the outer 
core, the position of the stalled shock. 

The loss of neutrinos at low and intermediate densities is important in 
determining the saturation ratio of the lepton to baryon ratio at the time 
of core bounce (the leptons determine the pressure in the core whereas the 
baryons, mainly nucleons, determine the mass of the homologous collapsing 
core). Brown et al [87] argued that the hydrodynamic collapse is nearly ho- 
mologous, i.e. the density structure of the collapsing core remains self-similar 
throughout the collapse until the time of bounce, this greatly simplifies the 
study of various processes during the collapse and quantities in the core can 
be calculated (largely semi-analytically) through the evolution of a "mean" 
"one-zone" symbolizing the core properties (see e.g. [88], [89]). 

The core of a massive star collapses under its own gravity when the pres- 
sure support from degenerate electrons is reduced through the capture of 
electrons in the stellar material. The electron capture on neutron rich heavy 
nuclei initially proceeds primarily through the allowed type (Al = 0) Gamow 
Teller transitions. As core density exceeds ~ 10 11 gm/cm 3 , the nuclei become 
become more and more massive and too neutron rich to allow e _ -capture to 
take place through allowed Gamow- Teller transitions from the ground state. 
This is because the allowed states for p to n transition within the nucleus 
are already filled by the neutrons (neutrons shell blocked) and the transi- 
tion strength for typical captures like 56 Fe — * 56 Mn used earlier (as in [90] ) 
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is no longer representative of typical nuclear e~-capture rates. It was shown 
[HI], [HI] that the dominant unique first forbidden transition strength was actu- 
ally negligible compared to the thermally unblocked strength under the typical 
core collapse conditions. Therefore, after neutron-shell blocking sets in, (when 
(A,Z) > 74 Ge) the sum rule for the Gamow Teller transition operator \Mqt\ 2 
decreases from a typical value of 2.5 [90], [88] to about 0.1. The e~-capture 
rate on a single nucleus X(A, Z) in the initial state i to the final state j is: 

A. i= l n 2^feM (5 6 ) 

J 13 

where fUj is related to |Mgt| 2 by ftij = 3.596 — log|McT| 2 for allowed 
Gamow- Teller type transitions (for free protons, log ft/. p . = 3.035). The fac- 
tor /i(T, ji e ,Qij) is the phase space factor for allowed transition, which is a 
function of the ambient temperature T, the Fermi-energy of the electron yU e 
and the Q- value for the reaction. The neutrino energy is E Ur — E e — Qij. 

The change in entropy during collapse controls the fraction of the dripped 
protons with respect to that of the heavy nuclei and this influences the overall 
neutrino spectrum received on earth as the spectrum of neutrinos generated 
by electron capture on protons are different from captures on heavy nuclei. 
The received neutrino spectrum depends not only upon the initial conditions 
from which the collapse started, but also on the details of the electron capture 
properties of the stellar matter. Properties of nuclei at finite temperatures and 
density during this phase of the collapse, where shell and pairing corrections 
are relevant were computed in 94J and utilized to evolve self-consistently 
with the electron capture physics and the consequent changes in nuclear and 
thermodynamic variables. 

The rate of generation of neutrinos per nucleon within energy band E v 
to E v + dE u after accounting for the relative abundance of free protons and 
nuclei, is: 

dY v {E v ) = d\ fp (E v )X p + d\ H (E v )(l -X n - X p )/A (57) 

here A represents the Atomic weight of the ensemble of nuclei present in 
the core, taken to be represented by a single "mean" nucleus as in [50]. The 
differential neutrino production rates for free protons and heavy nuclei are: 



dX log2 <G> El(E u + Q fpM )^/{E~TQ fp , H ^-(rn e c^) dE 

fp ' H (ft) fp .H (m e c 2 ) 5 (l + exp(E 1/ + Q fPj H -He)) 

(58) 

where the Coulomb correction factor < G > has been taken as « 2 for heavy 
nuclei and 1 for free protons. The Q-value is given as: Q = (/t+ 1.297+ Eqt) 
assuming that the strength is concentrated in a single state and here fx {— fi n — 
H P ) is the difference in the neutron and proton chemical potentials when free 
nuclcons coexist with a distribution of neutron rich nuclei in nuclear statistical 



Thermonuclear fusion and collapse and explosion in massive stars 



45 



equilibrium, and Eqt is the energy of the Gamow- Teller Resonance centroid. 
The centroids in fp-shell nuclei, found from experimental data from (n,p) 
reactions have been used for characterizing GT transitions in these nuclei [53] 
and are close to the value (3 MeV) used here. The GT centroid in the e-capture 
direction is itself a function of the ambient temperatures as a quasi-particle 
random phase approximation (QRPA) calculation shows [95]. At the relevant 
high densities, when neutron rich nuclei with A > 65, are abundant and the 
electron chemical potential is noticeably larger than typical nuclear Q-values, 
electron capture rates are mainly sensitive to the centroid and total strength 
of the GT+ distributions - these are reasonably well described within the 
random phase approximation [96] . We note that the temperature dependence 
of the nuclear symmetry energy can also affect the neutronization of the stellar 
core prior to neutrino trapping during gravitational collapse [97], [98] since 
ambient temperatures can reach upto several MeV. Not only the reaction Q- 
values but also the equation of state of bulk dense matter, the free nucleon 
abundances, the degree of dissociation into alpha-particles and the nuclear 
internal excitations are modified by changes in the symmetry energy. 

The Fermi-energy fi e = ll.l(pioy e ) 1/ ' 3 MeV. The difference in chemical 
potentials p,, and the relative fraction of free protons are obtained from a low 
density analytic equation of state similar to that in [90 with modifications 
noted in [88]. Shown in Fig. [15] is a typical "snapshot" spectrum of neutrinos 
from a 15 Mq star's core collapse within a narrow range of stellar core density 
around 10 11 gm cm~ 3 . Note the two separate peaks from captures on free 
protons and heavy nuclei which have non-thermal spectra. 

8.2 Number of neutrinos emitted and predictions of detections 

The number of neutrinos and their energy distributions that could have been 
detected by large underground neutrino experiments like the SuperKamioka 
experiment and the Sudbury Neutrino Observatory (SNO0 has been calcu- 
lated [5FJ. The spectrally integrated fluence of v e at a distance of 1 kpc, as 
Y e changes from 0.42 to 0.39 in a 1.4 M Q stellar core (of a 15 M Q star) is: 
F v = 4.2 x 10 1:L cm _2 . The energy of the infall neutrino burst up to this stage 
is: E Ve = 7.2 x 10 50 erg. The flux, direction and the spectra of the neutrinos 
could have been measured by the charge current dissociation of the deuterium 
nucleus (v e (d 7 pp)e~) in the ("classic") Sudbury Neutrino Observatory (SNO) 
[100 with a fiducial mass of 1 Kt of high purity D2O. It would have been 
also possible to detect v e and obtain spectral information by means of the 

20 The original SNO experiment with heavy water has finished taking data. Thus 
the heavy water based predictions if a nearby supernova took place are only 
indicative. However there are plans of an extended SNO (called SNO+) which 
will use 1 kTon of liquid scintillator which will greatly reduce the lower energy 
threshold of neutrinos and will remain sensitive to neutrinos from a nearby SN. 
Moreover, after an accident the photomultiplier tubes in SuperKamioka, the low 
energy thresholds originally designed for are no longer operative. 
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Fig. 15. "Snapshot" neutrino fluence (MeV 1 baryon ) in a density interval 
Ap 10 = 0.0002 around p 10 = 9.8668 gm cm" 3 for a 15 M star (D= 1 kpc) [9"9"] . 



neutrino-electron scattering reaction (y e + e~~ — > v e + e~) in SNO as well as in 
the light water detector Super Kamioka. Apart from the reactions mentioned 
above which can perform spectroscopy (i.e. measure the incoming neutrino en- 
ergy), the neutral current dissociation of d by the reaction v e {d 1 pn)v e could 
have obtained the the total neutrino flux (of all flavors). The number of v e 
events which could have been detected in the SNO detector through neutrino- 
induced c.c. reaction on the target d nuclei is given by: n Vi! = F v a cc {e Ve )N,i 
where Nd (= 6.02 x 10 31 ) is the total number of target nuclei present in the 
lKt detector. The charge current and neutral current cross-sections {a cc {e Vc ) 
and er nc (e„ x ) respectively), have been computed for the v-d process by Bahcall 
et al. [101] and accurate fits to these cross-sections between 5 to 40 MeV are 
given as [102] : a{ — on{t v — e t h,i) 2 ' 3 where i = cc and nc, a cc = 1.7 x 10~ 44 
cm 2 , a nc = 0.85 x 10" 44 cm 2 , e thjCC = 2.2 MeV and e th ,nc = 1-44 MeV. 

For the H2O based Cerenkov detector (Super-Kamioka) the (v e ,e~) scat- 
tering events would be the main source of v e spectral information since the 
corresponding energy thresholds for charge current and neutral current inter- 
actions for ordinary water are much higher. (During the collapse stage, the 
neutrino flux is almost entirely in neutrinos of the electron type; anti-neutrinos 
of various kinds, as well as neutrinos of the mu or tau type are generated in co- 
pious numbers only in the hot post core bounce phase). The relevant (f e , e~) 
scattering cross-section is [TUg] : er e = (l/2)(4G 2 m 2 ?l 2 /7rc 2 )(7/12)(e l/ /TO e c 2 ). 
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The number of detections by Super-Kamioka (mass 32 kt) and SNO for a 
supernova explosion 1 kpc away, for several possible scenarios of stellar core 
collapse are reported in Table 2 [53] ■ The 15 M Q star collapse is initiated 
from thermodynamic conditions as in [55J (Y e j = 0.420, Si/ks = 1.00, T, = 
0.7178, pio = 0.37), while the 25 M Q star's single zone initial conditions are 
similarly derived from the data reported in [104 and an expression for the 
core averaged entropy (Y el = 0.423, S l /k B = 1.14, T* = 0.6756, p w = 0.15). 

Table 2. Pre-trapping neutrino detections predicted in SNO (heavy water) and Su- 
per Kamioka |99| with hardness ratios up to pio = 24.16 for indicated heavy nuclear 
e-capture matrix elements for 15 Mq Fuller (1982) and 25 Mg WWF presupernova 
stars. Note the caveat in footnote 20. 



Star Mass 


\M G t\ 2 


tcollapse 


Pre-trappin 


g Variables 


No. Detected 


Hardness Ratio ^ | 






(ms) 


Yef 


Sf/ka 


SNO 


S-K 


SNO S-K 


15 M 


1.2/0.1 


120 


0.3909 


1.0021 


82 


394 


0.2786 0.8540 




2.5/0.1 


120 


0.3896 


1.0085 


66 


344 


0.2876 0.9537 


25 Mq 


1.2/0.1 


190 


0.3828 


1.1080 


120 


566 


0.2878 0.8319 




2.5/0.1 


190 


0.3813 


1.1204 


99 


499 


0.2916 0.9190 



' The hardness ratio denotes the number of neutrino events in the 5 MeV < E Ve < 12 
MeV and 12 MeV < E Ve < 25 MeV bands. 



Neutrino spectroscopy of the final state of a star would be possible provided 
that the event occurs at a relatively close distance. Although, a priori, these 
may be rare events, there have been a number of historical Supernovae, as 
well as detected radio pulsars within a distance of about 2 kpc. There are a 
number of star forming regions nearby (such as the Orion complex - about 
440 pc away), which are sites of potential supernova progenitor. A detection 
by underground neutrino experiments would constrain features of theoretical 
calculations of both collapse and explosion era of type II Supernovae as well as 
shed light on the characteristics of the stellar core. The 19 neutrinos detected 
from SN1987A were most likely to have been emitted during the post-bounce 
phase as their total flucncc during the proto neutron star cooling phase (at 
~ 10 58 ) is much larger than that during the collapse phase(~ 10 56 ). 

9 Detected neutrinos from SN 1987A and future 
neutrino watch 

Gravitational collapse of the core of the massive star under its own grav- 
ity leads to a supernova explosion. These are extremely energetic explosions 
where the observable energy in the kinetic energy of the exploded debris and 
electromagnetic radiation add up to several times 10 51 erg. The actual energy 
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scale is typically 3 x 10 53 erg or higher, but most of this is radiated away 
neutrinos. Although the full understanding of the process of explosion in a 
gravitational collapse leading to a supernova has not been achieved despite 
several decades of theoretical and computational work, a watershed in the field 
was achieved observationally when a supernova exploded close by in a satellite 
galaxy of our own, namely SN1987A in the Large Magellanic Cloud (LMC). 
A few neutrinos were detected from this supernova [105], [106] . which were 
the first detections of the astrophysical neutrinos from outside of our solar 
system. By using the energetics of the neutrinos, their approximate spectral 
distribution, the distance to the LMC it was possible to show that the overall 
energy of the explosion was indeed Et ~ 2 - 3 x 10 53 erg. 

In addition, the duration of the neutrino burst was of the order of a few 
seconds as opposed to a few milliseconds, signifying that the observed neu- 
trinos had to diffuse out of the dense and opaque stellar matter as predicted 
theoretically, instead of directly streaming out of the core. The spectral char- 
acteristics also indicated that the object that is radiating the neutrinos was 
indeed very compact, roughly of the same dimensions as that of a protoneu- 
tron star. Thus SN1987A provided the observational confirmation of the broad 
aspects of the theoretical investigation of stellar collapse and explosion. For a 
review of the understanding of the astrophysics of SN1987A, see [107] . 

Physicists are now gearing up to detect not only another supernova in our 
own galaxy, but by hoping to build very large neutrino detectors, they aim 
to detect supernova neutrinos from the local group of galaxies ( |108j . |109j ). 
As neutrinos from the supernova travel directly out from the core, they arrive 
a few hours ahead of the light flash from the exploding star, since light can 
only be radiated from the surface of the star, and it takes the supernova 
shock launched at the deep core several hours to reach the surface. In the 
case of SN1987A, this time delay allowed the estimation of the size of the star 
that exploded. Thus some advance warning of the supernova in the optical 
band can be gotten from a "neutrino watch" . Physicists have now connected a 
worldwide array of neutrino detectors through the Internet (SN Early Warning 
System or SNEWS^j]) which will notify astronomers to turn their optical, UV 
and other telescopes to the right direction when they find a burst of neutrinos 
characteristic of a supernova explosion. 

10 What X-ray spectroscopy reveals about 
nucleosynthesis in SNe and SNRs 

X-rays from a supernova explosion arise from the interaction of the super- 
sonic ejecta with the circumstellar medium (CSM). The CSM typically con- 
sists of a slow-moving wind. When the ejecta collides with the CSM, it creates 
two shocks: a high-temperature, low-density, forward-shock ploughing through 



See the site: 



http://hep.bu.edu/~snnet / 
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the CSM (known as blast-wave shock) and a low-temperature, high-density, 
reverse-shock moving into the expanding ejecta. Initially the X-rays come from 
the forward-shocked shell dominated by continuum radiation, but after a few 
days X-rays arise also from the reverse-shock, which can have substantial line 
emission, thus providing nucleosynthetic fingerprints of the ejecta. The tem- 
poral evolution of the X-ray luminosity of a supernova can yield information 
on the density distribution in the outer parts of the exploding star (p oc ?'~™, 
here n can be in the range 7 — 12, typically n ~ 10 for a Blue Supergiant 
(BSG) and n ~ 12 for a Red Supergiant (RSG)- see |110j ). These studies are 
therefore of interest to stellar structure and evolution. 

To date thirty-six supernovae have been detected in the X-ray banda^l. 
Among these the most extensively studied is SN 1987A because it was so 
bright. Another extensively studied object, of the supernova remnant (SNR) 
kind is Cassiopeia A. A SNR is just an older supernova after it has picked 
up enough material from the surrounding medium, which slows it down. The 
debris of the explosion, or the ejecta, radiates when it is reheated after an 
initial cooling, having been hit by a inward propagating reverse shock. 

10.1 Supernova Remnant Cassiopeia A 

Cassiopeia A (Cas A for short and also known as 3C461 or Gill. 7-2.1) is 
the second youngest SNR in our Milky Way galaxy. It was believed to be the 
product of a SN explosion in ~ 1672 1112] only about 3.4 kpc away |113j . The 
British astronomer John Flamsteed may have recorded it as a sixth magni- 
tude star in 1680, but it may have faded rapidly after explosion which could 
have acted against its widespread reportage. It is a shell type SNR which was 
rediscovered by radio astronomers Ryle and Smith [114j and is the brightest 
extra-solar radio object. Cas A is an oxygen-rich SNR with heavy element 
distribution typical of a core-collapse SN. Its closeness, young age and high 
brightness across the whole electromagnetic spectrum have underscored its 
importance for studying supernova^!. A recent optical spectrum of the orig- 
inal supernova near maximum brightness, obtained from scattered light echo 
more than three centuries after the direct light of the explosion was received 
on Earth, showed it to be a type lib SN |116) 24 l . somewhat like the well studied 



See S. Immler's X-ray supernova page at 

http:/ /lheawww. gsfc.nasa.gov/users/immler/supernovae_list. html and [111] , 
Here we discuss thermal emission from Cas A, which directly connects to nu- 
cleosynthetic products. Many shell-type SNRs, including Cas A, also show non- 
thermal synchrotron X-rays - see [115] . 

Krause et al note that even with the overall lack of hydrogen emission in most 
knots and the nitrogen enrichment in the remnant which is widely interpreted 
as signatures of the collapse of a Wolf-Rayet star, i.e. a type lb SN |117j . Cas 
A cannot be classified as arising out of a type lb SN, since Cas A light echo 
spectrum does not match well the spectrum of the prototype SN 2005bf. 
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prototype SN 1993J0 [123 . The estimated mass loss rate of 2 x 1(T 5 M yr" 1 
and a wind velocity of 10 km s _1 consistent with the hydrodynamical state 
of the Cas A remnant |124j is also similar to what has been interpreted from 
radio and X-ray observations of SN 1993 J [125], [125] . 

SNR Cassiopeia A has a bright clumpy emission ring with a diameter of 
about 3'. This is associated with the SNR ejecta, whereas at an outer diameter 
of about 5', there is fainter emission in a filamentary ring which is due to the 
forward ("blast-wave") shock in the circumstellar medium. The SNR has a 
nearly circular appearance, though the ejecta shows a bipolar structure, and 
jets (mainly in the NE, i.e. upper left of the image and a fainter extension in 
the SW direction, called the "counterjet") that were formed by the explosion 
(detected both in the optical [127] and in the Si line X-ray emission [129J, 
[128] ). The remnant has been extensively observed in the optical and more 
recently with a million second exposure of the Chandra X-ray Observatory 
|128] (see also XMM- Newton [130] ) and includes high resolution grating X- 
ray spectroscopy with Chandra [131j as also the Spitzer Space Telescope using 
the Infrared Array Camera 132J and the IR Spectrograph 133 . 

Many of Cas A's optically bright knots have been identified as shocked 
ejecta which is still visible due to its young age. The optical radiation emitting 
regions have been classified into two groups: the so called fast moving knots 
(or 'FMK's with speeds 4000 km s _1 < v < 15000 km s" 1 ) and the slow 
moving quasi-stationery floculi (or 'QSF's with speeds v < 300 km s _1 ). The 
FMKs believed to comprise of SN ejecta, have H-deficient emission (lacks in 
Ha) which is dominated by forbidden lines of O and S emission [134| while 
the QSF emission is rich in N, and is believed to originate from the gas lost 
from the star before its explosion (circumstellar envelope) which is now hit 
by the blast wave shock [131j . 

Cas A shows the nucleosynthesis products of both hydrostatic and explo- 
sive nucleosynthesis. C and O are produced in He burning, Ne and Mg first 
appears through C burning, and O and Al are added with Ne burning [135] . 
Spectral lines of heavier elements: Si, S, Ar, Ca etc are seen. These are pro- 
duced in O burning, with alpha reactions on Mg also contributing, - they are 
from zones where explosive O-burning and incomplete explosive Si-burning 
occurs. Fe group elements are produced in the "Si burning" chain of reactions 
and result from complete and incomplete Si-burning. Much of the layered 

25 The collapse of a red supergiant [118] in a binary system [119] with the larger 
star with a mass on the main sequence of ~ 10 to 20 led to SN 1993J and 
could explain its light curve and other characteristics [120] . [121] , Photometry 
and spectroscopy at the position of the fading SN with Hubble Space Telescope 
and Keck Telescope, a decade after its explosion showed the signature of a massive 
star, - the companion to the progenitor that exploded. While the binary system 
initially consisted of e.g. a 15Mq star and 14M© star in an orbit of 5.8 yr, at the 
time of the explosion, the primary had a mass of 5.4Mq (with a helium exhausted 
core of 5.IM0) and a secondary which gained mass in transfer and which we still 
see today, ended up with 22M© [122] , 
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nucleosynthetic structure has been preserved during the explosion process in 
Cas A, with layers of N-, S-, O-rich ejecta seen beyond the outer shock |134j . 
However, in the optical bands, "mixed emission knots" showing N and S lines 
suggest that clumps of high speed S ejecta have penetrated through outer 
N-rich layers [134]. X-ray line data show that for example, in the southeast 
region of the SNR, Fe is farther out and is moving faster than the Si/O region. 
Hughes et al [136 consider the high surface brightness knots enriched in Si 
and S are consistent with nucleosynthetic products of explosive O burning, 
whereas elsewhere, the more diffuse, lower surface brightness features with 
enhanced Fe could be the result of explosive Si burning. 
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Fig. 16. Spectrum of the entire supernova remnant Cassiopeia A taken by Chandra 
(AXAF CCD Imaging Spectrometer (ACIS)) on January 30-31, 2000 against pulse 
height (lower x-axis) in the spectral channel and photon energy (upper x-axis) using 
a gain value of 4.8 eV per channel (Fig. reproduced by permission of the AAS and 
courtesy of U. Hwang |129] b 
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The "non-dispersed'^ X-ray CCD spectrum for the entire Cas A remnant 
is shown in Fig. [TH] A 50 ks exposure in January 2000 collected approximately 
16 million photons during the observation of the SNR. It shows prominent 
lines of Si, S, Ar, and Ca, in their Helium-like ionic state and undergoing 
(n = 2 — > n = 1) transitions as also the L and K transitions of Fe. 

Maps were constructed of the entire remnant in the spectral line energies 
of several specific elements. The Si and S X-ray line maps and the optical 

26 This means the X-ray beam did not pass through any dispersive element like a 
grating, but the spectrum is obtained by measuring the energy of the incident 
X-ray by the intrinsic spectral sensitivity of the CCD. 
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Fig. 17. Top: Three-color image of Cas A with red = Si Hea (1.782.0 keV), blue 
= Fe K (6.526.95 keV), and green = 4.26.4 keV continuum. The remnant is roughly 
5' across. Bottom left: Overexposed broadband image showing faint features. The 
spectral regions are indicated (top left of this box: northeast jet; bottom left of 
this box: Fe-rich region; lines at bottom right point to two southwest jet filaments). 
Bottom right: On the same scale, the ratio image of the Si Hea (1.782.0 keV) and 
1.31.6 keV (Mg Hea, Fe L), without subtraction of the continuum contribution. The 
image highlights the jet and counterjet traced by Si emission, although features at 
the lowest intensity levels are uncertain. (Figure reproduced by permission of the 
A AS and courtesy of U. Hwang |128p . 

maps resemble each other and this demonstrated that the X-ray and optically 
emitting ejecta are largely spatially coexistent. In fact the X-ray Si and S line 
emission, expected to include a significant contribution from the ejecta, shows 
asymmetric Doppler shifts corresponding to bulk velocities of ~ 2000 km s _1 
which are comparable to those of the optical ejecta knots [129j . Si, S, Ar, and 
Ca have similar general morphologies, but examination of their images in their 
line energies also show significant variations. At the positions of the brightest 
knots in the Si image at the inner boundary of the shell to the northeast and 
southeast the corresponding Ar and Ca X-ray images are much weaker in the 
remnant. Hughes ct al [136 using the Chandra first light observation of Cas A 
(~ 5000s on 1999 August 20) showed that the Fe-rich ejecta lie outside the S- 
rich material, and claimed that this is due to extensive, energetic bulk motions 
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which caused a spatial inversion of a part of the supernova core. However, on 
the basis of Spitzer Space Telescope imaging and spectral data Rudnick and 
collaborators [137] . |133| claim that there are two roughly spherical shocks, 
(the blast wave and the reverse shocks) but "the wide spatial variations in 
composition seen in infrared, optical and X-ray studies are not due to local 
differences in ionization age or temperature, but instead reflect very specific 
asymmetries in the geometry of the underlying explosion. These asymmetries 
are ones in which the nucleosynthetic layers have remained mostly intact along 
each radial direction, but the velocity profiles along different radial directions 
vary over a range of approximately five" (see also |132j ). In some directions, 
only the upper C-burning layers have been probed by the reverse shock, while 
in other directions, deeper O- and Si-burning layers have reached the reverse 
shock. 
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Fig. 18. Composite Spitzer image of Cas A in [Ar III] (blue), [S IV] (green), and 
[Ne II] (red), showing two pronounced neon-rich crescent-shaped regions to the N 
and S. The size of the image is 320" with north to the top and east to the left. Shown 
are the kinematic center of the remnant at the magenta circle and X-ray localization 
of the remnant's compact object (red square, 7" to the south). Green lines indicate 
the la range of the "kick vector" direction of the compact object from the ejecta's 
expansion center: 169° ±8.4° [112] . The two regions of enhanced neon abundance lie 
very close to this projected direction. Several much smaller neon-enhanced regions 
lie in the West along the X-ray jet direction (Fig reproduced by permission of AAS 
and courtesy of J. Smith [133] ). 
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The gas infrared lines of Ar, Ne, 0, Si, S and Fe are seen at different 
locations, along with higher ionization states of the same elements visible in 
the optical and X-ray lines. For example, strong [Ar II] regions in the IR match 
very closely the places where helium- and hydrogen-like ionization states of 
Si and S are seen in the X-ray and [S III] seen in the optical [134) . These 
are the same regions with multiple temperatures and ionization states but 
show products of oxygen burning instead of carbon burning. Other regions, 
for example the crescent-shaped green feature slightly east of south in the 
Spitzer IRAC line and continuum images depicting gas in Fig 2 of [132] show 
strong [Ne II] emission in the IR (see also the red regions in the image Fig [TBI 
reproduced from [133) ). increasing [O III] emission in the optical and a gap in 
the silicon dominated X-ray emission, are interpreted as the locations where 
carbon-burning layers are presently encountering the reverse shock. Different 
layers containing various types of nucleosynthetic products seems to show the 
presence of characteristic types of dust, e.g. the deep layers contain silicates 
while the upper layers contain dust dominated by AI2O3 and carbon grains. 
They also find evidence for circumstellar dust heated by the blast wave shock 

Ml- 

Smith et al [133j identify IR line emission from ejecta materials in the 
interior, prior to their encounter with the reverse shock, as well as from the 
post-shock bright ring. There is a dramatic increase in the electron density (by 
a factor > 100 to ~ 10 4 cm~ 3 ) as well as a concomitant change in the ionization 
state of the ejecta as it encounters the reverse shock. There is a clear layering of 
ionization state, from low ionization species in the interior, e.g. [Si II] (8.2 eV), 
higher energies on the IR bright ring [S IV] (34.8 eV) where optical emission 
is also seen, and very high energies traced by X-ray line emission from H- 
like and He-like K-alpha resonance lines of Si XIII (0.5 keV) and Si XIV (2.4 
keV), extending beyond the IR-bright rim. In addition, they find two compact, 
crescent shaped clumps with high neon abundance (mentioned above) which 
are arranged symmetrically around the central neutron star, and the crescent 
regions are closely aligned with the kick direction of the neutron star from the 
remnant's expansion center. These regions contain a huge amount of ionized 
neon (dominated by Ne + and Ne ++ , and excluding neutral neon), almost 
1.8 x 10~ 4 M Q , flowing outwards 20 degrees from the plane of the sky at 
roughly —5500 km s _1 in the south and +4200 km s _1 in the north, while 
the entire SNR may contain an ionized neon mass ~ 8.6 x 10 -4 M Q . Smith 
et al comment that the apparent macroscopic elemental mixing mentioned in 
|136j . may actually arise from different compositional layers of ejecta passing 
through the reverse shock at present along different directions. 

X-ray grating spectra of Cassiopeia A and SN 1987A 

In the optical bands, thousands of individual knots have been observed to yield 
kinematic information, whereas the X-rays probe more dynamically important 
information since a much larger fraction of the ejecta mass is probed by X-rays 
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(~ 4M as compared to < 0.1M Q in the optical). The Chandra High Energy 
Grating Spectroscopy (HETGS) [131j provided for the first time the high 
spatial and spectral resolution X-ray map of Cas A, comparable to what is 
achievable in the optical. Because of Cas A's bright emission lines and narrow 
bright filaments and small bright clumps that stand out against the diffuse, 
continuum emission, this grating X-ray spectroscopy was possible and yielded 
rich information about the kinematic and plasma properties of the emitting 
knots. The HETGS has two grating arms with different dispersion directions, 
with the medium energy range (MEG) covering the range 0.4 — 5.0 keV and 
spectral FWHM of 0.023 A, while the high energy range (HEG) covers 0.9 — 
10.0 keV and spectral FWHM 0.012 A. The grating spectra were obtained by 
Lazendic et al on some 17 different positions on the Cas A image, mainly near 
the reverse shock regions which are X-ray bright. 

As mentioned already, Cas A image in the HETGS both in the non- 
dispersed (zeroth-order) as well as dispersed images in different energy bands 
contain H- and He-like ionization states of O+Nc+Fe (0.65-1.2 keV), Mg 
(1.25-1.55 keV), Si (1.72-2.25 keV), S (2.28-2.93 keV), Ar (2.96-3.20 keV), Ca 
(3.75-4.0 keV). In addition images in the Li- and Be-like states of Fe K lines 
(6.3-6.85 keV) are also obtained. A key advantage of high resolution grating 
spectroscopy over the low resolution "non-dispersive" CCD spectroscopy is 
that the former can be used to resolve individual lines and thereby interpret 
conditions in the radiating plasma in various parts of the SNR. The He-like 
ions of Si and S and other elements present in Cas A are the dominant ion 
species for each element over a wide range of temperatures and they emit 
strong K-shell lines of these ionic stages. For typical plasma densities in the 
SNR, the He-like triplet of Si and S lines shows strong forbidden (f) and res- 
onance lines (r) and a comparatively weaker intercombination line (i) (see 
section 4.5 of [138 for an illustration of these transitions in O-ions). The G- 
ratio = (f+i)/r is, for example, a good diagnostic of the electron temperature 
|139] . Thus the measured line ratios and abundance ratios of H- and He-like 
ions of the same element (say Si) were used |131j to measure the electron 
temperature kT e and the ionization time scale r = n e t (n e being the elec- 
tron density and t being for example the time since the region was hit by 
the shock). The results of [131] show that for most of the selected regions, 
resolved spectroscopy of Si He-like triplet lines and Si H-like lines yields tem- 
peratures around ~ 1 keV, consistent with reverse shocked ejecta. However, 
two regions, (designated R8 and R10) had significantly higher temperatures 
~ 4 keV, and could be part of the circumstellar material. The dominant el- 
ement in Cas A seems to be Oxygen, and the electron density in the X-ray 
line emitting regions seem to vary between 20 to 200 cm -3 , which is more 
typical of the unshocked interior of the IR line-emitting plasma, than that 
at the shock front as shown in 1 133 (see above). Lazendic et al had derived 
unambiguous Doppler shifts for their selected 17 regions. While the SE region 
of the SNR show mostly blue shifted velocities reaching upto —2500 km s^ 1 
(compare the IR line velocities in the south —5500 km s^ 1 ), the NW side of 
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the SNR had extreme red shifted values up to +4000 km s _1 (compare IR 
line velocity in the north: +4200 km s _1 ). 

Deep high resolution X-ray spectroscopy of SN 1987A was undertaken 
with Chandra HETG, 20 years after its explosion |140j . The impact of the SN 
ejecta with the circumstellar medium dominates the observed luminosity, the 
X-ray luminosity having brightened by a factor of 25 from that in 1999 and is 
presently increasing at the rate of 40 % per year [141] . An expanding elliptical 
ring is seen in the X-ray image whose brightness distribution seems to correlate 
with the rapidly brightening optical hot spots on the inner circumstellar ring 
observed with Hubble Space Telescope. The HETG spectrum for a total live 
time of 355 ksec shows H-like and He-like lines of Si, Mg, and Ne, and O VIII 
lines and bright Fe XVII lines. Fig. [T9l taken from Dewey et al |140] shows 
the resolved Si and Mg triplets of (r, i, and f lines); the data have similarity 
with model G-ratios (/ + i)/r mentioned earlier. Since the dispersed X-ray 
spectrum is a convolution of the spatial structure of the X-ray image and the 
motion of the X-ray emitting gas, Chandra's sub-arc second resolution was 
important to resolve the circumstellar ring of dimensions 1.2" x 1.6". 




6.6 6.7 6.8 9.1 9.2 9.3 

Wavelength [Angstrom] Wavelength [Angstrom] 



Fig. 19. Si (left panel) and Mg (right panel) triplets (resonance, intercombination 
and forbidden lines) resolved by Chandra High Energy transmission Grating Spec- 
troscopy. The data are shown by the solid black histogram and an arbitrarily scaled 
point-source version of the two-shock model (no spatial- velocity effects) is shown in 
red to suggest the relative similarity of the data and model G-ratios [(f+i) /r] . (Fig 
reproduced by permission of the A AS and courtesy of D. Dewey [140| h 

The global fit to the HETGS data with a two shock model yielded element 
abundances and absorption column densities which are consistent within the 
same 90% confidence limits that were derived from LETGS by |142j . The lower 
temperature shock seems to give the same kTi ow for all the LETG data sets (in 
2004 and 2007) and the HETG observation, whereas the kThigh values seem 
to show a general evolution towards lower values among the grating datasets 
and similar trend seen in the ACIS monitoring of SN 1987A. 
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The HETG data implies a relatively low bulk radial velocity of the shocked 
gas in the ring compared to model expectation from a plane parallel strong 
shock entering stationary gas and the temperature range inferred from the 
spectral modeling of the emitting gas. Since the X-ray image is correlated 
with the optical hot spots, it appears that the blast wave ahead of the SN 
ejecta is overtaking dense clumps of circumstellar gas in the hot spots. The 
blast wave may be encountering the dense clumps either at normal incidence 
or at oblique incidence. In the former case, the reflected shock would give rise 
to gas that has been shocked twice and having nearly stationary bilk velocity 
but further elevated temperature. With more and more X-ray emission coming 
from gas behind reflected shocks, there would be an increase of the fraction of 
higher temperature X-ray emission, consistent with what is seen. If the blast 
wave encounters dense clumps at an oblique incidence, the shocked gas will 
have significant velocity component parallel to the shock surface, and a signif- 
icantly fast transmitted component emitting X-rays can result. The Doppler 
broadening seen in the line profiles of the Chandra data seems to suggest that 
both transmitted and reflected shocks encountering the circumstellar ring at 
normal and oblique angles may be at play [140] . 

10.2 Live radioactive decays in Cas A, SN 1987A 

Gamma-ray, X-ray, optical and infrared line spectroscopy of SNe and SNRs 
have been used to observe nucleosynthesis and the abundance distribution of 
elements freshly synthesized (both radioactive isotopes and stable decay prod- 
ucts) and to extract dynamical information about the explosion. In particular 
the radioactive isotopes provide unique tracers of nucleosynthetic processes 
(what, where and how much) and its related dynamics. The best examples 
are the observations of gamma-ray lines in supernovae, but X-ray line spec- 
troscopy of decay products of radioactive nuclei have also been attempted, 
and specific elements in numerous SNRs and a few SNe identified. 

After a few days, the main energy input to the SN ejecta comes from 
radioactive decay. For SN 1987A, at first, the important isotopes are 56 Ni 
(t^ = 6.1d) followed by 56 Co (r 1/2 = 78.8(2). Beyond - 1100 days, 57 Co is 
more important and at > 2000c£, the dominant role is played by 44 Ti. Using 
bolometric and broad-band UBVRI light curves [143] estimated the masses 
of the three most important radio-isotopes in SN 1987A to be 56 A r i(0.069 ± 
0.003M Q ), 57 A^(0.003M Q ) and 44 Ti(l ± (0.5 - 2) x 10- 4 M Q ). 

Among the radioactive isotopes accessible to gamma-ray astronomy, 44 Ti 
is a key isotope for the investigation of the inner regions of core collapse 
SNe and their young SNRs. It has a half-life of 58.9 ± 0.3 yr [144] and a 
decay scheme: Ai Ti ^ 44 Sc ^ 44 Ca. The discovery [145] of the 1157 keV 
gamma-ray line emission^ 7 ! from Cas A with the Compton Gamma Ray Ob- 
servatory was the first direct proof of synthesis of this short lived, freshly 

27 = 2 + — » + transition in 44 Ca reached from 44 Sc following electron capture. 
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made radio-isotope in SNe. Renaud et al [146] using the INTEGRAL space- 
craft recently reported the detection of both 67.9 and 78.4 keV gamma-ray 
lines of 44 <Sc in Cas A. There was a clear separation of the two lines due to 
an improved detection of the hard X-ray continuum up to 100 keV. The line 
flux of (2.5 ± 0.3) x 10~ 5 cm~ 2 s _1 leads to a tightly constrained 44 Ti mass 
of (1.6lo 3) x 10 -4 -Mq). This is actually high compared to the predictions of 
the standard models [135] , [147j , or their improved versions [148] , [149] . Since 
the production of 44 Ti is sensitive to the explosion energy and asymmetries 
and Cas A is known to be asymmetric and energetic (2 x 10 51 erg instead 
of the standard 1 x 10 51 erg), this could be a factor in its apparent overpro- 
duction compared to models. At the same time, the recent revision of the 
40 Ca(a, 7) 44 Ti reaction rate |150j has led to an increase of 4i Ti production 
by a factor of ~ 2. 

Live radioactive isotopes freshly synthesized in the explosion were detected 
from SN 1987A by directly detecting 1238 keV and 847 keV gamma-ray lines 
[15 LJ , which provided detailed diagnostics of the nuclear burning conditions 
and the explosion dynamics. Another radio-nuclide 57 Co, which is a decay 
product of 57 Ni made in the supernova explosion, was also detected directly 
by measuring gamma-ray lines [152J. The calculated gamma-ray light curves 
for 847 keV and 1238 keV and other lines |153j could be made consistent with 
the Solar Maximum Mission (SMM) measurement only if the 56 Co was mixed 
up to a mass coordinate of M r ~ 13.5M Q . The calculated fluxes were still 
slightly smaller than the observed fluxes at t ~ 200d, but this discrepancy 
could be removed by taking account of the effect of dumpiness in the ejecta 
which is more significant in the earlier phases. Also, the observed flux ratio 
between the two gamma-ray lines at 847 keV and 1238 keV is close to unity 
at early stages because of the smaller cross section for the 1238 keV lines 
than for 847 keV. As the column depth of the overlying matter decreased 
due to dilution with time in an expanding envelope, the observed flux ratio 
approached the expected value of 0.68 based on branching ratios [153] . 

X-ray fluorescence spectrometry technique has been used by archaeologists 
to determine the major and trace elements in ceramics, glass etc. Here, a sam- 
ple is irradiated with X-rays and the wavelengths of the released (fluorescent) 
X-rays are measured. Since different elements have characteristic wavelengths, 
they can be identified and their concentrations estimated from the intensity of 
the released X-rays. Trace elements analysis may, for example, help in iden- 
tifying the (geographical) source of a material. A similar method has been 
proposed [154j and attempted for SN 1987A 155 . It utilizes electron capture 
decays of freshly synthesized radio-nuclei. Following the electron capture, the 
K-shell vacancies are in most cases filled in by downward transitions from 
other bound shells. As the fluorescence X-ray yields are well known, mea- 
suring the X-ray line fluxes can then estimate the number of freshly synthe- 
sized nuclei. For example, attempts were made to detect the 5.9 keV Ka line 
from the stable nucleus 55 Mn in the Chandra X-ray spectrum of SN 1987A, 
which is due to the decay of radioactive 55 Fe. Both 55 i 7, e(r 1 / 2 = 2.7 yr) and 
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55 Co(r 1 / 2 = 18 hr) are produced in explosive, incomplete Si burning as well 
as in normal freezeout of nuclear statistical equilibrium, in the inner ejecta of 
core collapse supernovae. However, no evidence of the 5.9 keV line emission 
from Mn could be found in 400 ks of Chandra ACIS data and the upper limit 
to the mean flux was < 3 x 10~ 7 cm _2 s _1 . Rauscher et al |148j calculated the 
ejected mass of A — 55 radioactive nuclei to be 7.7 x 10~ 4 Mq for 20M© mod- 
els of which most was 55 Co. If only about half this mass of 55 Fe were ejected, 
the reduced flux would be consistent with the observed upper limit. On the 
other hand the, even if the total mass inside were as much as 1 x 10~ 3 M Q , but 
the 55 Fe abundance was zero outside the radial velocity shells at 1500 kms -1 , 
the line flux would be still consistent with data, as at late times the emerging 
flux depends sensitively on the presence of 55 Fe in the outer zones. 

10.3 Other X-ray supernovae 

Apart from SN 1987A, which occurred close by and was therefore easy to 
detect, there are supernovae with relatively high intrinsic X-ray luminosity. 
SN 1993J and SN 1995N are two of them and they are at the high end of 
the X-ray luminosity |156j . [157] . which makes it suitable for study in the 
X-ray wave bands even at late stages |159j . Nymark et al [158] have shown 
that SNe with strong X-ray emission are likely to have radiative shocks and 
that in these shocks a large range of temperatures contribute to the spectrum. 
The cooling of the shock also affects the hydrodynamics of the flow structure 
and the volume of the emitting gas and thus the total luminosity from the 
interaction region. A shock is radiative depending on whether the cooling time 
scale t coo i of the shock is short compared to its expansion time t. Since the 
cooling function (from electron bremsstrahlung or free-free scattering and X- 
ray line emission from bound-bound transitions of ionized atomic species) is a 
function of the shock temperature this condition translates to a condition on 
the shock temperature for radiative vs. adiabatic shock. The boundary of the 
two regimes depends upon the shock velocity V s and the ejecta density profile 
(given by the index rj in p oc v^ 71 ) and the separation in turn depends upon 
the chemical composition of the ejecta (see Fig. [20l reproduced from [160J ) . 
For a given shock velocity (and a composition and mass loss parameter of the 
progenitor stellar wind and the ejecta velocity scale), the shock is likely to be 
radiative for steeper density gradients. 

As already mentioned, the spectroscopic type of SN 1993J from a type II- 
like SN at early times due to detected Balmer lines soon weakened and it came 
to resemble a type lb SN, leading to a reclassification as a type lib SN |161j . 
This transition is best explained if the progenitor star prior to its explosion, 
had lost its hydrogen envelope due to interaction with a binary companion 
[120] . [121], [162], [163], [122]. If the H-envelope was indeed already thin, 
then the reverse shock would transit the H-rich region quickly, and soon the 
emission may be dominated by material dredged up from the nuclear processed 
parts of the star in the interior. In the 4H47 model of Shigeyama et al [164] , by 
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Fig. 20. The border between radiative and adiabatic regions of the shocks in 
supernovae as a function of shock velocity and ejecta density profile index 77 (p oc 
v -vy -pjjg boundaries depend upon the composition of the ejecta through which the 
shock may be propagating as also the mass loss parameter of the progenitor stellar 
wind and the ejecta velocity scale (Fig. courtesy T. Nymark [160| V 

2600 days after explosion, the hydrogen-rich envelope had ~ 0.47Af Q . Inside 
this envelope is the outer part of the He-zone which is enriched in N (about 
~ 0.3M Q ), while its inner part is C-rich (has ~ 1.5M Q ). If the circumstellar 
medium or the ejecta is clumpy, they may be hit by the reverse shock obliquely 
and the shock in the clump will be slower and the shock temperature lower. 
The presence of both hard X-rays and optical emission from the shocked CSM 
shows the importance of a clumpy CSM, for example in SN 1987A where a 
range of shock speeds are necessary to explain the observed spectrum |142j . 
SN 1993J is not as bright as SN 1987A or Cas A due to its larger distance, 
and both Chandra and XMM-Newton data is comparatively of a moderate 
signal to noise, at the low (CCD) resolution. The XMM spectrum of SN 1993J 
is dominated by the Fe L emission at 0.7 - 1.0 keV and is blended with strong 
Ne IX-X emission. Above 1.0 keV, emission lines which may be present are 
from Mg XI-XII, Si XIII-XIV, as well as S XVI. Nymark et al find that the 
spectrum of SN 1993J is fit best by a combination of an adiabatic shock 
{kT rev = 2.1 keV) propagating through a zone with CNO burning products, 
with He dominant and N and Ne being the most abundant metals, and a 
second radiative shock (kT rev =1.0 keV). The latter could have been caused 
by instabilities at the reverse shock front or by a clumpy CSM. The XMM 
spectrum of SN 1993 J and the model fits used [160j are shown in Fig|2"Tl 

SN 1995N is a type Iln supernova, a type that shows unusual optical char- 
acteristics and spans a very broad range of photometric properties such as 
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Fig. 21. XMM spectrum of SN 1993 J (shown by data-points with crosses from 
April 2001). The data is fitted with an adiabatic shock at kT rev = 2.1 keV (long 
dashed in blue) and a radiative shock at kT rev — 1.0 keV (short dashed in red) for 
a He/N dominated composition from model sl5s7b [163] . (Fig. courtesy T. Nymark 
H60]). 



decline rates at late times |165| . It is likely that these differences are related 
to their progenitor's structure, mass, composition as well as the composition 
and the density profile of the CSM [166] . These supernovae show the pres- 
ence of strong, narrow Balmer line emission on top of the broader emission 
lines in their early spectra. The narrow emission lines may originate in the 
dense and ionized circumstellar (CS) gas [1671 1168] . The presence of strong 
Ha emission line, the high bolometric luminosity and the broad Ha emission 
base powered by the interaction of the supernova shock with the CSM, all 
point towards a very dense circumstellar environment |169j . When the stellar 
core collapses and explodes, the supernova lights up the slow-moving gas into 
narrow emission lines leading to the type Iln supernova classification (n for 
narrow emission line). This interaction of the supernova shock with the dense 
CS gas is indicated by strong radio and X-ray emission detected from several 
type Iln supernovae and in particular SN 1995N. 

Chandra ACIS observations of SN 1995N show a best-fit line energy at 
1.02 keV which is ascribed to Ne X [170] . There is also the possibility of 
another line with best-fit line energy at 0.85 keV with identification as Ne IX. 
The ionization potentials of all ionized Neon species upto Ne VIII are less 
than 240 eV, whereas those of Ne IX and Ne X are above 1195 eV. Hence at 
temperatures found by the Chandra observations, the predominant species of 
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Neon are expected to be Ne X and Ne XI. The mass of Neon was estimated 
to be about 5 x 10~ 3 Af Q — 1 x 10~ 2 M Q . This is consistent with Neon being 
in the Helium layer for a 15M Q star where it was co-synthesized with C, N, 
and He. For other stellar masses and other composition zones, the required 
Neon mass is much larger than observed. Therefore, Neon in the Helium core 
of the 15M© star is the probable site where it was co-synthesized with C, N. 

Conclusions 

Much insight about how stars burn nuclear fuel, evolve and ultimately explode 
has been obtained by calculations and computer simulations. Theoretical de- 
velopments have been augmented by crucial astronomical observations over 
decades in all bands of electromagnetic radiation and through other signal 
channels like neutrinos. Most important have been the wider availability of 
powerful telescopes both from the ground and space platforms and state of 
the art computers. Crucial inputs to the field of nuclear astrophysics are also 
coming from laboratory experiments involving radioactive ion beams (RIB) 
and intense beams of energetic nucleons and nuclei. Short lived nuclei can 
only be studied close to their sites of formation in the laboratory before they 
decay. Such facilities will further advance the future of nuclear astrophysics. 
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